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ABSTRACT 


The  law  of  conservation  of  momentum  is  interpreted  for  an 

i 

electromagnetic  system.  Using  Maxwell’s  stress  tensor  and  the 
concept  of  electromagnetic  angular  momentum  flow  an  expression  is 
obtained  to  calculate  the  instantaneous  electromagnetic  torque  on 
a  body  in  free  space  with  a  view  to  apply  it  to  the  theory  of 
electrical  machines. 

The  electromagnetic  field  problem  of  the  double  cylindrical 
machines  is  reduced  to  finding  the  transient  solution  of  a  vector 
wave  equation.  The  method  of  the  Laplace  transform  is  used  to 
develop  a  general  solution  of  the  equation.  The  solution  is  exact 
and  may  include  any  number  of  space  harmonics  and  time  functions 
associated  with  the  driving  frequencies  and  the  characteristic  roots 
of  the  electromagnetic  wave  system. 

The  machine  volt-ampere  equations  are  obtained  by  integrating 
the  electric  field  intensities  and  the  torque  equation  is  deduced 
from  the  momentum  concept.  Quasi-static  field  quantities  and  the 
resulting  machine  equations  of  motion  are  also  derived  from  the  low 
argument  approximations  of  the  Bessel  functions  involved  in  the 
exact  equations. 

An  exact  wave  theory  analysis  of  an  induction  machine  is  devel¬ 
oped  in  detail  to  study  the  production  of  torque  in  terms  of  a  trans¬ 
fer  of  electromagnetic  angular  momentum  across  the  air  gap.  The 
simpler  case  of  the  machine  under  steady  state  conditions  is  treated 
first  and  subsequently  a  solution  is  given  for  a  particular  type  of 


- 


i 


. 


transient.  It  is  found  that  the  quasi-static  approximation  provides 
completely  satisfactory  quantitative  information  in  the  steady  state, 
but  the  exact  treatment  is  necessary  to  give  a  proper  account  of  the 
torque  at  each  instant  during  the  course  of  the  transient. 
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CHAPTER  I 


INTRODUCTION 

1.1  Historical  Development  of  the  Concept  of  Electromagnetic 

Momentum 

Newton’s  corpuscular  theory  of  light  led  its  adherents  in  the 
early  eighteenth  century  to  associate  momentum  with  radiant  energy. 
Many  investigators1  in  that  century  directed  powerful  beams  of  light 
on  delicately  suspended  bodies  and  tried  to  detect  a  pressure  due  to 
the  impulse  of  the  corpuscles.  Motions  were  observed  but  the  experi¬ 
ments  were  inconclusive  due  to  thermal  effects  which  were  not  properly 
understood  at  that  time.  From  a  theoretical  study  Euler  inferred 
the  existence  of  light  pressure  in  1748.  He  further  accounted  for 
the  tails  of  comets  by  supposing  that  solar  rays,  impinging  on  the 
atmosphere  of  a  comet  drive  off  the  more  ^subtle  of  its  particles. 

In  the  late  nineteenth  century,  J.C.  Maxwell  gave  his  celebrated 
electromagnetic  theory  which  regards  light  and  radiant  heat  as 
phenomena  of  electromagnetic  wave  notion.  On  the  basis  of  his  theory 
he0  deduced  in  1873  that  in  a  medium  in  which  waves  are  propagated 
there  is  a  pressure  in  the  direction  normal  to  the  waves  and  numer¬ 
ically  equal  to  the  energy  density.  Thus,  if  a  beam  of  radiation  is 
incident  perpendicularly  on  a  perfectly  reflecting  metallic  sheet  the 
radiation  pressure  is  equal  to  the  density  of  energy  in  the  region 
immediately  outside  the  metal.  In  1876  Adolfo  Bartoli-*-  argued  that 
if  radiant  energy  is  transported  from  a  cold  body  to  a  hot  one  by 


*7'  -% 


\  '  ■  :  .V  -  " 


-1 


, 


: 


'  y-/: .  j'.'  'V-  r-fi.-> ■  •  :  ■  *  V  •  •"• 


-•  .  ' '1  •  -Si 

..  •  '  :  '  '  ■  •{'»  '  ; 


'  *  -  ■. 


i 

■ 

; 

. .  .  V 

■ 

\ 

MA 

••It 

- 

Jt 

.. 

1>  ■  ' 

J 

'4 

J 

... 

’’ 

i 

. 

• 

. 

. 

/  '  .  ■ 

:  -  ;  ■'•: 

.. 

,  -  • 

■ 

• 

-  .  ''  ;  -  ■ 

' 

• 

•  - 

1,  . 

4 

•  r  ^ 

■  ' 

a  -  ‘v 

■ 

.  . . .  * 

. 

• 

• 

f„  •  • 

* 

*  •- 

a  ir 

’  . 


; 


'  V.  ■■  .  ' 

\  *■ 

.  .»  ■  '  &  m  ■  -■  ■  ■ 


■■  ' :  .:.r ..  • 


...  ?  ■:  ■  •  '  ■ 


•• .  i 

"  • 


2 


means  of  a  moving  mirror,  the  second  law  of  thermodynamics  would  be 
violated  unless  a  pressure  were  exerted  on  the  mirror  by  the  light. 

In  1900,  Planck’s  quantum  theory1*  further  reinforced  the  associa¬ 
tion  of  momentum  with  radiation. 

The  first  reliable  experimental  verifications  of  Maxwell’s 

theoretical  result  on  radiation  pressure  were  done  by  Lebedew^ 

in  1899  and  E.F.  Nichols  and  G.F.  Hull3  in  1901.  N.  Carrara  and 

P.  Lombardini7 ,  in  1949,  detected  radiation  pressure  of  centimeter 

waves  by  using  a  725A-type  magnetron  and  a  rectangular  wave  guide 

•  •  ft  9 

system.  In  1952,  Cullen  and  his  associates  ’  developed  an  appara¬ 
tus  for  the  absolute  measurement  of  microwave  power  in  terms  of 
radiation  pressure.  Their  instrument  consisted  essentially  of  a 
small  vane,  suspended  in  rectangular  wave  guide  carrying  an  mode. 

If  linear  momentum  could  be  associated  with  a  plane  polarized 
wave,  it  was  a  small  step  to  associate  angular  momentum  with  a  circu¬ 
larly  polarized  wave.  This  was  first  done  by  Sadowsky13  in  1899, 
and  later  by  J.H.  Poynting11  in  1909  and  Epstein13  in  1914.  With 
such  an  angular  momentum  they  also  predicted  an  associated  torque  on 
a  material  body.  Dr.  Poynting  compared  the  angular  momentum  with 
that  of  a  revolving  shaft. 

In  1935,  R.A.  Beth13  passed  circularly  polarized  light  through 
a  doubly  refracting  plate  which  changed  the  state  of  polarization 
and  he  was  able  to  detect  a  torque,  thus  making  the  first  experimen¬ 
tal  verification  of  the  existence  of  electromagnetic  angular  momentum. 
N.  Carrara14  reported  experimental  works  at  microwave  frequency  in 
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1949.  Recently,  at  the  University  of  Alberta,  F.S.  Chute^  measured 

within  experimental  error  the  reaction  torque  on  an  antenna  that 

radiated  a  net  angular  momentum  at  radio  frequency. 

16 

Historically,  aether  played  an  important  role  in  man’s  imagina¬ 
tion  in  developing  different  philosophies,  and  its  model~  varied  from 

age  to  age.  The  momentum  concept  associated  with  electromagnetic 

17 

theory  was  first  mathematically  developed  by  Sir  J.J.  Thomson  ,  a 
follower  of  Maxwell,  in  the  last  decade  of  the  nineteenth  century. 

He  adopted  the  view  that  aether  was  the  vehicle  of  momentum  in  appar¬ 
ently  free  space  and  in  an  electromagnetic  field  problem  the  aetheral 
momentum  must  be  taken  into  account  for  Newton's  third  law  of  motion 
and  the  principle  of  conservation  of  momentum  to  be  valid. 

The  hypothesis  that  aether  is  a  store  house  of  momentum  was 

17  18 

further  developed  by  Poincare  5  ,  who  in  1900,  referring  to  Thomson's 

2 

expression  that  in  free  aether  the  electromagnetic  momentum  is  (1/c  ) 
times  the  Poynting  flux  of  energy,  suggested  that  electromagnetic 

O 

energy  might  possess  mass  density  equal  to  (1/c  )  times  the  energy 
density . 

The  failure  of  the  Michelson-Morley  experiment  led  to  the 
abandonment  of  the  idea  of  a  material  aether  and  resulted  in  the 
ultimate  development  of  the  special  theory  of  relativity  by  Einstein. 

An  essential  feature  of  this  theory  is  the  association  of  mass  with 
energy  and  the  enormous  success  of  the  theory  is  probably  the  strong¬ 
est  argument  in  favour  of  the  concept  of  electromagnetic  angular  momentum. 
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1.2  Survey  of  the  Field  Theory  of  Electrical  Machines 

Machine  engineers  have  attempted  from  time  to  time  to  adapt 
Maxwell’s  equations  to  study  different  aspects  of  electrical  machines. 
Modern  generalized  machine  analysis  includes  the  formulation  of 
machine  differential  equations  in  certain  reference  system  in  terms 
of  the  lumped  parameters  of  circuit  theory  and  the  solution  of  these 
equations  for  specified  load  and  terminal  conditions.  The  quasi¬ 
static  solution  of  Maxwell’s  equation  gives  the  approximate  electro¬ 
magnetic  fields  in  different  parts  of  the  machine.  The  magnetic 
fields  thus  obtained  are  utilized  to  define  inductance  parameters 
through  the  concept  of  energy  storages  or  flux  linkages.  The 
differential  equations  are  then  formed  by  using  basic  force  laws 
or  Euler- Lagrange  equation  of  classical  mechanics. 

For  specified  conditions  the  machine  equations  may  be  treated 
by  various  topographical  methods-^  including  an  equivalent  circuit 
approach  for  easy  cases  or  by  the  more  generalized  tensorial  meth- 
ods^5^.  These  methods  may  collectively  be  called  the  ’’system 
methods”  and  most  of  the  time  there  is  an  interplay  between  them. 

The  quasi-static  field  components  are  sometimes  used  in  connection 
with  Poyntings  vector  to  calculate  the  power  transfer  through  the 
air  gap  of  the  machine. 

OO 

Historically,  perhaps  Creedy  was  the  first  to  discuss  wave 
phenomena  in  the  dynamo-electric  machine  in  1917.  Hague ^  published 
a  series  of  papers  between  1917  and  1926,  where  he  discussed  electro¬ 
magnetic  fields  and  forces  in  various  regions  of  a  machine. 
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Ollendorff  was  the  first  to  use  vector  potentials  in  machine  theory. 

The  application  of  Poynting's  vector  to  energy  calculations  in 
low  frequency  power  apparatus  and  systems  was  first  given  by  Slepiarh^ 
in  1919.  He  tried  to  give  a  picture  of  the  way  the  insulating  medium 
supporting  electromagnetic  waves  plays  a  vital  role  in  leading  the 

9  c 

electrical  power  to  where  it  is  consumed.  In  1942,  Slepian^  again, 
in  a  very  ingenious  derivation  of  the  energy  flow  theorem  gave  nine 
alternative  forms  of  the  power  flow  vector.  Alger  and  Oney^7  s  28 
studied  Poynting’s  energy  flow  in  induction  machines  and  Hawthorne f 
did  the  same  in  synchronous  machines.  Alger  and  Oney  tried  to  visua¬ 
lize  the  machine  as  an  array  of  magnetized  air  spaces  in  which  intense 

and  turbulent  energy  concentrations  exist. 

30 

In  1954,  Mishkin  obtained  approximate  solutions  of  Maxwell’s 
equations  for  an  idealized  model  of  an  induction  motor.  His  model 
replaced  the  toothed  stator  and  rotor  by  continuous,  but  inhomogeneous 
and  anisotropic  regions,  having  averaged  properties.  Cullen  and 
Barton f  analyzed  the  induction  machine  in  1958,  by  making  use  of 
the  wave  impedance  concept  and  employing  a  transmission  line  analogue. 
The  sleeve-rotor  machine  was  analyzed  by  Guilford ^  in  1962.  Besides 
solving  Laplace's  equation  in  the  air  gap,  he  solved  the  diffusion 
equation  for  the  conducting  sleeve. 

Professor  Saunders'^?  34  contribution  to  the  field  theory  of 
electrical  machines  is  quite  significant.  Using  probable  practical 
conductor  distributions  in  a  double  cylindrical  structure  he  generated 
surface  harmonic  current  sheets  and  obtained  quasi-static  solutions 
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of  Maxwell's  equations  for  tlie  air  gap  fields.  In  his  field  equations 
there  were  provisions  to  include  time  harmonics  in  the  surface  curr¬ 
ents.  A  Lagrangian  formulation  was  then  used  to  find  the  equations 

35 

of  motion.  Recently,  Nasar  published  a  review  paper  on  the  electro¬ 
magnetic  theory  of  electrical  machines ,  giving  a  brief  digest  of  exist¬ 
ing  information  in  a  unified  form  by  classifying  diverse  methods  of 
analysis.  A  large  bibliography  on  the  subject  is  also  provided. 

1. 3  Background  of  the  Electromagnetic  Angular  Momentum  Concept 

Applied  to  Induction  Machine  Theory 

31 

In  1958,  Cullen  and  Barton  mentioned  that  the  theory  of  the 
squirrel  cage  induction  machine  is  closely  related  to  the  phenomenon 
of  radiation  pressure.  They  tried  to  draw  an  analogy  between  the 
force  acting  on  the  rotor  of  an  induction  motor  and  the  radiation 

stress  of  a  surface  wave,  obliquely  incident  on  an  absorbing  surface. 

36  37 

In  two  recent  papers ,  Professor  Barlow  5  refers  to  the  induc¬ 
tion  motor  as  a  device  in  which  torque  is  produced  by  the  absorption 
of  electromagnetic  angular  momentum.  The  inference  is  that  electro¬ 
magnetic  angular  momentum  is  generated  by  the  stator  winding  and  is 
propagated  across  the  air  gap  to  the  rotor,  in  much  the  same  way  as 
angular  momentum  is  propagated  in  a  waveguide  supporting  a  circularly 
polarized  wave. 

The  above  three  references  mention  the  mathematical  relation 
P  =  Too,  where  P  is  the  power,  T  is  the  torque  and  go  is  the  angular 
velocity  of  the  field  components.  At  this  point  it  is  appropriate  to 
note  that  this  relation  is  very  well  known  to  induction  machine  engin- 
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eers ;  but  they  derived  it  without  ever  associating  momentum  properties 
with  the  electromagnetic  field.  The  concept  of  electromagnetic  angular 
momentum  in  relation  to  induction  machine  torque  production  has  not  yet 
been  developed  in  the  mathematical  sense  and  it  is  a  problem  to  be  dis¬ 
cussed  in  this  thesis. 

1.4  Conservation  of  Energy  and  Momentum  in  Electromagnetic  Systems  - 
Poynting,s  Equation  and  Maxwell's  Stress  Tensor 
Two  of  the  fundamental  laws  of  classical  physics  are  the  conserv¬ 
ation  of  energy  and  momentum.  Since  energy  and  momentum  of  material 
bodies  can  be  affected  by  an  interaction  with  an  electromagnetic  field, 
it  is  convenient  and  proper  to  associate  an  electromagnetic  energy  and 
an  electromagnetic  momentum  with  the  field  in  such  a  way  as  to  pre¬ 
serve  the  conservation  laws. 

In  the  case  of  energy,  a  few  simple  manipulations  of  Maxwell’s 
equations  lead  to  the  famous  Poynting’s  equation,  namely 


. . .1.1 


°S 


The  second  term  on  the  right-hand  side  of  equation  (1.1)  is 
interpreted  as  the  rate  of  increase  of  electric  and  magnetic  energy 
within  the  volume;  the  first  term  includes  the  power  dissipated  as 
Joule  heat  and  any  other  electromagnetic  energy  conversion.  From 
the  conservation  of  energy  the  left-hand  side  of  equation  (1.1)  then 
represent  the  inward  flow  of  electromagnetic  energy  across  the  sur¬ 


face  bounding  the  volume. 

A  similar  profit  and  loss  account  can  be  made  for  momentum.  Con- 
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sider  the  electromagnetic  volume  force  density  given  by 


38 


f=  ptE  +  JtxB  ...1.2 

•  •  • 

where  is  the  sum  of  free  and  polarization  charges  and  J^_ 

includes  free,  polarization  and  amperian  current  densities.  The 
electromagnetic  force  density  given  by  equation  (1.2)  may  be  equated 
to  the  time  rate  of  change  of  the  mechanical  momentum  density  g^, 


that  is: 


dgm  _  p  E  +  ?xB 
dt  11  t 


... 1. 3 


Using  Maxwell’s  equations  it  may  be  shown  that 

P  E  +  J  *B  =  %E  V.E  +  +  €0Bx|| 

L  L  c)Z 

since  £oBx||  =  £0Ex||  -  &.(  €0  ExB) 


. . .1.4 


and 


one  may  write 


=  -  £0ExV*E  -  (  P-0EaB) 

di 


tv.  % 

^ o 


=  0 


f  +  A(4oe,S)  =  c, 


EV.E  -  EX7/E 


+A 


^  V-  B  - 


. . .  1 .  S 


Introduce  the  symmetric  tensor  M  of  rank  two,  known  as  Maxwell's 
stress  tensor,  defined  in  dyadic  form  by 


M  =  P^E  +  -±-  tg  -  k  £0E.E  +  ±4.g)U 

/*o  z  Mq 


... 1. 6 
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where  U  is  the  fundamental  tensor  whose  components  are  ,  the 
Kronecker  delta. 

In  appendix  (1)  it  is  shown  that  the  right-hand  side  of 
equation  (1.5)  is  equal  to  the  divergence  of  M,  where 


div  M  =  (M^)  SJ 
Hence  it  follows  that 


1,2,3. 


.1.7 


div  M 


f  + 


at 


60e*b) 
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Taking  the  volume  integral  of  both  sides  of  (1.8)  and  using  the 
divergence  theorem  of  Gauss  applied  to  tensor  of  rank  two,  one  may  write 


M.ndS  = 


^5 


(RE  +  J  xB ) dV  + 

v  t  ^  °  « 


(  £0E*B)dV 


1.9 


Let  the  surface  S  now  recede  to  infinity  or  at  least  far  enough 
to  yield  a  negligible  contribution  so  that  the  left-hand  side  of 
equation  (1.9)  is  zero.  Combining  this  with  equation  (1.2)  one 
obtains 


Therefore , 
The  quantity 


dt 


»(gm  + 


^m 


V 
^m  = 


£0ExB)dV  =  0 
^SxBMV  -  Constant 
£oEXB 
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is  dimensionally  a  momentum  per  unit  volume  and  may  be  referred  to  as 
the  volume  density  of  electromagnetic  linear  momentum.  In  an  electro¬ 
magnetic  system,  the  sum  of  mechanical  and  electromagnetic  momentum 
is  conserved  —  the  mechanical  momentum  is  not  necessarily  conserved. 

There  is  a  marked  similarity  between  equations  (1.1)  and  (1.9). 
The  first  integral  on  the  R.H.S.  of  (1.9)  is  the  total  mechanical 
force  which  by  Newton's  second  law  is  the  rate  at  which  mechanical 
momentum  is  produced.  Since  the  integrand  is  expressed  in  terms  of 
electrical  quantities  it  may  be  interpreted  as  the  rate  at  which 
electromagnetic  momentum  is  converted  to  another-  form  of  momentum. 

The  second  integral  on  the  R.H.S.  gives  the  rate  of  increase  of 
electromagnetic  momentum  within  the  volume.  Consequently,  by  the 
law  of  conservation  of  momentum  the  L.H.S,  gives  the  inward  flux  of 
electromagnetic  momentum  through  the  bounding  surface ,  S . 

39 

If  Maxwell's  equations  are  written  in  symmetrical  form  by 
defining  a  magnetic  current  density  Jj?  =  (/^ M)  and.  a  magnetic 

charge  density  p*  =  -div^jM,  then  the  force 


f  =  pE  +  JXyU0H  +  p*H  -  J*  X  £0E 


. . . 1. 2A 


where 


P  =  Pf  -  V- P 


Then 


J  =  Jf 


M  =  €0EE  +  -  y(^E.E  +y%H.S)U 


1.6A 


and 


1.11A 


The  difference  in  the  two  definitions  of  momentum  density  given 
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by  (1.11)  and  (1.11A)  is  due  to  the  fact  that  the  mechanical  force 
density  acting  inside  a  material  medium  has  been  given  by  different 
expressions.  There  would  appear  to  be  no  experimental  procedure  by 
which  this  problem  can  be  resolved.  It  is  important  to  note  that 
in  free  space  the  two  approaches  are  identical.  In  this  thesis  only 
the  overall  forces  on  material  bodies  (rotor  and  stator)  are  con¬ 
sidered,  and  not  the  forces  existing  within  the  bodies.  As  a  result 
field  quantities  existing  in  the  air  gap  only  are  considered  in  dis¬ 
cussing  the  transfer  of  momentum. 

Since  the  linear  momentum  and  angular  momentum  are  conserved 
independently,  in  general,  the  electromagnetic  field  will  have  an 
angular  momentum  with  respect  to  certain  point  or  axis.  Take  the 
moment  of  equation  (1.8)  such  that 


...1.12 


where  r  is  the  radius  vector  from  the  arbitrary  origin  to  the  field 
point.  Integration  over  the  volume  concerned  yields 


1.13 


Again  r  x  div  M  dV  has  the  important  property  that  the  integral 
over  any  volume  can  be  transformed  into  an  integral  over  the  surface 
S  enclosing  the  volume  such  that 
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|  f’xdiv  MdV 

ly 


rx(M.n)dS 


. . .1.14 


Hence 


frx(M.ft)dS  =  r*x?dV  +  ^  M)dV 

Jv  V  ...1.15 


Equation  (1.15)  expresses  the  conservation  of  angular  momentum  in 
an  electromagnetic  system.  The  left-hand  side  is  the  total  inward 
flux  of  electromagnetic  angular  momentum  through  the  surface.  'The 
first  term  on  the  R.H.S.  gives  the  mechanical  torque  produced,  that 
is,  the -rate. of  conversion  of  electromagnetic  angular  momentum  to 
mechanical  angular  momentum.  The  second  term  may  be  interpreted  as 
the  rate  at  which  .the  electromagnetic  angular  momentum  stored  within 
the  volume  is  changing.  As  in  the  case  of  linear  momentum, 

geu  = 

is  taken-  as- the-  electromagnetic  angular  momentum  density.  From 
equation  ( 1 . 11A) . 

geu  = 

The  flux  density  of  the  angular  momentum  of  a  circularly  polar¬ 
ized  wave  about  the  axis  of  propagation  may  be  shown  to  be  P/to,  where 
P  is  the  power  flow . per  -unit  area  and  to  is  the  angular  frequency. 

This  relation  may  also  .be  derived  from  a  simple  quantum  argument. 

Each  photon  comprising  the  beam  has  an  angular  momentum  +  Ti/2tt,  where 
h  is  . Planck' S:  constant..  If  w  is  the  energy  density,  the  total  num¬ 
ber  of  photons  .per  unit  volume  is  n  =  w/hf  =  P/chf ,  where  f  is  the 
frequency  of  radiation.  If  the  spins  of  all  the  photons  are  aligned 
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in  the  positive  sense, 

£>em 


the  total  angular  momentum  per  unit  volume  is 

P _ h_  _  P_ 

chf  2tt  ooc 


This  momentum  travels  at  the  velocity  of  light  and  hence  the  flux 
density  of  angular  momentum  crossing  unit  area  normal  to  the  direction 

p  P 

of  propagation  is  equal  to  If  the  waves  are  absorbed  by  a 

OOC  03 

screen  disposed  perpendicularly  to  the  direction  of  propagation,  the 

p 

torque  per  unit  area  would  be  -. 

1. 5  Electromagnetic  Torque  on  a  Body  in  Free  Space 

In  an  electromagnetic  system  if  a  body  is  surrounded  by  free 
space  (air) ,  then 


rx(M.n)dS 


gives  uniquely  the  total  inward  flow  of  electromagnetic  angular  mom¬ 
entum.  This  has  been  discussed  in  the  previous  section. 

Since  for  sinusoidal  steady  state  the  time  averages  of 


9  r  -*  -*  8 

—  ?x(£  ExB)dV  and  — 

3t  L  °  3t 


->  -* 


rxyo€Q(ExH)dV 

V 


are  zero,  the  mechanical  torque  of  electromagnetic  origin  on  the  body 
can  be  obtained  uniquely  from 


T  =  frx(M.n)dS  ...1.16 

em  Jg 

(S  in  free  space) 

For  the  transient  state,  the  full  equation 


P 

JV 


rxfdV  = 


rx(M.n)dS  -  — 

r>  3t 


SeodV 
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.  . .1.17 


has  to  be  used. 


'  -h;  -V  . 


. 

'  •  t  • 


'  f* 


- 


■  :  < 

- 


■a 


...  r 

'  V 


...  r 

. 


i 

4 


[f 


i  .  - 

v,.  .  .•  'i 


«#!•  +■ 


“  ‘  '  '  '  ' 


'  '  ■  \ 
\  \  i. 

< ;  .  •  *  >  ’  •  - 


r 


■i  ■' 


*■ 


. 


' 


n 


14 


If  for  physical  reasons  the  electric  and  magnetic  field  intensiities 
inside  the  body  are  vanishingly  small,  the  momentum  storage  inside  it 
is  negligible.  For  such  a  body  the  mechanical  torque  for  both  steady 
and  transient  state  is  given  by 

T  =  f?x(M.ft)dS  -  (?  x  (M)dV  ...1.17A 

em  Js  dt  b-V 

0 

where  V  =  the  volume  enclosed  by  S 
VQ  =  volume  of  the  body 

This  equation  will  be  applicable  in  an  induction  machine  (in  fact 
in  any  machine)  with  the  rotor  having  iron  of  very  high  permeability. 
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CHAPTER  2 


THE  MACHINE  MODEL  AND  THE 
BOUNDARY  VALUE  PROBLEM 

2.1  Introduction 

Field  theory  analysis  of  electrical  machines  involves  solutions 
of  appropriate  differential  equations  derived  from  Maxwell's  equa¬ 
tions  for  specified  conditions.  In  the  mathematical  sense  it  is  a 
boundary  value  problem  with  due  consideration  of  the  initial  condi¬ 
tions  for  transient  analysis.  Normal  machine  geometry  is  so 
complicated  that  exact  analysis  of  the  electromagnetic  field  problem 
is  highly  complex  and  sometimes  impracticable.  This  requires  the 
reduction  of  the  actual  machine  into  a  simpler  mathematical  model 
with  some  appropriate  and  reasonable  assumptions. 

2 . 2  Cylindrical  Model 

In  the  simplest  and  quite  common  machine  model,  the  inner 
surface  of  the  stator  and  the  outer  surface  of  the  rotor  are  concen¬ 
tric  cylinders,  separated  by  a  small  air  gap.  On  these  surfaces 
are  slots,  parallel  to  or  slightly  skewed,  in  which  lie  the  current 
carrying  conductors.  The  distribution  of  the  conductors  around  the 
periphery  of  the  machine  is  such  that  it  forms  a  periodic  function. 
The  conductors  belonging  to  the  different  phases  of  the  same  mag¬ 
netic  member  are  similar,  but  displaced  from  each  other  by  suitable 
angles.  The  stator  and  rotor  may  have  similar  winding  arrangements 
(wound  rotor  machines).  Quite  often  squirrel  cage  windings  are  used 
on  the  rotor.  Sometimes  a  continuous  conducting  element,  rather  than 
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a  network  of  individual  bars  is  used  on  the  rotor.  Such  a  rotor  may 
be  formed  by  plating,  which  will  provide  a  conducting  layer  on  the 
surface.  This  machine  may  be  called  a  sleeve  rotor  machine. 

2 . 3  Reference  Frame 

Because  of  the  cylindrical  symmetry  it  is  suitable  to  use 
cylindrical  co-ordinates  (r,0,z)  for  the  analysis  of  such  a  machine. 
Furthermore  it  is  convenient  to  use  reference  frames  fixed  to  each 
magnetic  member  such  that  the  stator  reference  system  is  stationary 
with  respect  to  stator  and  rotor  reference  system  is  stationary 
with  respect  to  rotor,  but  rotates  in  space  with  respect  to  the 
fixed  reference  frame  of  the  stator.  The  velocities  involved  are 
very  small  compared  to  that  of  light.  This  makes  time  and  length 
measure  in  the  two  systems  approximately  the  same  (no  relativistic 
time  and  length  difference) . 


Fig.  2.1  -  Schematic  Diagram  of  a  Wound  Rotor  Machine 
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Fig.  (2.1)  shows  the  schematic  diagram  of  a  two  phase  wound  rotor 
machine.  Only  the  centre  turns  of  each  phase  are  shown.  is  the 
space  angle  between  the  axes  of  stator  phase  (a)  and  rotor  phase 
(a).  When  the  rotor  rotates  <f>  is  a  function  of  time.  Therefore 
the  effect  of  rotation  must  be  taken  into  account  when  quantities 
of  one  frame  are  referred  to  another. 

2.4  Reduction  to  Boundary  Value  Problem 

Consider  Fig.  (2.2)  where  regions  (1),  (2)  and  (3)  represent 
respectively  air,  a  conducting  layer  with  continuous  current  distri¬ 
bution  and  iron. 


Region  (3)  -  Iron 
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Fig.  2.2  -  Three  -  Region  Boundary  Surfaces 
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The  sides  MN  and  PQ,  each  of  length  Al ,  of  the  contour  C  lie  in 
either  face  of  the  conducting  layer  and  the  sides  NP  and  QM  pene¬ 
trating  the  layer  are  of  length  As  each.  Integrate 

VxS  =  Jf  + 


over  the  surface  bounded  by  the  contour  and  apply  Stokes’  theorem 
to  obtain 


From  the  above  equation  it  may  be  shown  that  in  the  limit  as 
As-K),  i.e.  region  (2)  vanishes  to  a  surface 


T,  *(.%!,)=  Lim  (  ...2.2 

AS+O 

Because  D  and  its  derivatives  are  bounded,  the  first  term  on  the 
R.H.S.  of  equation  (2.2)  vanishes.  The  quantity 


1.  =  Lim 

2  AS-K> 


amp/m 


. . .2.3 


known  as  the  linear  current  density,  can  be  other  than  zero  on  the 
assumption  that  as  Thus 


t  Y  f  U  —  i 


...2.4 


For  a  distributed  winding  let  the  contour  C  with  sides  NP  and 
QM  in  the  tooth  regions  enclose  the  <th  slot.  Also  let  there  be  N 

lv 

conductors  in  this  slot  each  carrying  a  current  i  .  The  total 
current  flowing  in  the  slot  is  then  i  N  .  In  view  of  the  fact  that 
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y  -*30,  the  magnetic  field  intensity  along  NP,  PQ  and  QM  is  zero.  Hence 

the  line  integral  of  H  taken  around  the  slot  is  H  Al,  which  from 

Maxwell’s  equations  is  equal  to  the  current  in  the  slot.  Thus,  ft^A£  = 

N  i  .  If  the  number  of  turns  per  unit  lengfh  is  N.  =  N  /A£,  it 
k  k  ik  k 

follows  that  H_  =  N  n  i.  .  Clearly  i_  =  N-.,  i,  acts  as  a  linear 
1  lk  k  J  1  lk  k 

surface  current  density. 

For  the  distributed  winding,  the  air  gap  surface  may  now  be 

viewed  as  consisting  of  a  series  of  a  finite  number  of  discontinuous 

Linear  current  density  pulses.  Due  to  the  periodic  distribution  of 

these  pulses ,  it  is  possible  to  expand  them  into  Fourier  space  harm- 
33 

onics  .  For  the  continuous  winding  the  linear  current  density 

would  be  a  continuous  periodic  function.  In  fact  any  practical 

winding  including  the  effects  of  skewing  may  be  handled  in  this 
33 


manner  . 

The  electromagnetic  system  of  the  machine  is  thus  reduced  to  a 
boundary  value  problem  with  air  gap  surfaces  represented  by  Fourier 


series  of  linear  current  density.  Inside  iron  H-*0,  and  in  crossing 
the  surface  the  tangential  component  of  the  magnetic  field  changes 
discontinuously  by  an  amount  equal  to  the  linear  current  density. 

The  other  boundary  condition  may  be  shown  to  be  that  the  normal  com¬ 
ponent  of  the  magnetic  flux  density  must  be  continuous  across  any 

•  •  • 

surface  of  discontinuity.  This  follows  from  the  fact  that  V.B  =  0. 

2 . 5  The  Two-Phase  Model  to  be  Analysed  and  Summary  of 

Basic  Assumptions 

(a)  A  two  pole  machine  with  two  stator  windings  placed  in 
space  quadrature  will  be  considered.  It  will  be  taken 
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that  the  conductors  forming  a  pole  are  similar  to 
those  forming  the  other  pole  so  that  the  linear  current 
density  will  be  given  by  the  odd  harmonics , 


•  s 

1lna 


\  ,Ts  .s  .  ns 
/  N:  i  sin  n0 

_ ,n  a 


...2.5 


(b) 


(c) 


•  s 

1lnb 


NS  i^  sin  n  ( 0S— :tt/ 2 )  ...2.6 

-ji  b 


n=l,3 , . ; 


•  s  s 

where  i  and  i,  are  currents  in  stator  phase  (a)  and  (b) 

0.  D 

S 

respectively.  It  may  be  shown  that  is  a  constant 
corresponding  to  each  harmonic,  and  it  takes  care  of  the 
space  distribution  of  the  conductors. 

The  rotor  will  have  either  two  phase  windings  similar  to 
those  of  the  stator  or  a  sleeve  rotor  winding. 

The  machine  axial  length  is  considerably  larger  than  its 
radius.  (This  will  make  the  effect  of  the  currents  flow¬ 


ing  in  the  end  wires  negligible . ) 

(d)  The  relative  permeability  of  iron  is  taken  to  be  infin¬ 
itely  larger  than  that  of  air. 
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CHAPTER  3 


VECTOR  WAVE  EQUATION 


AND  BOUNDARY  CONDITIONS 


3.1  Introduction 

Exact  solution  of  Maxwell's  equations  are  required  to  establish 
the  concept  of  electromagnetic  angular  momentum  flow  in  an  induction 
machine.  The  solutions  may  be  attempted  by  defining  suitable  pot¬ 
ential  functions.  Because  of  the  assumed  linearity  of  the  system 
superposition  will  be  applicable.  It  has  been  stated  before  that 
each  reference  system  is  stationary  with  respect  to  the  phase 
winding  producing  the  field.  Hence,  the  mathematical  form  of  the 
field  equations  and  their  solutions  for  each  phase  will  be  similar. 

3. 2  The  Air  Gap  Vector  Potential 

In  appendix  (2)  Maxwell's  equations  are  solved  in  terms  of 
suitably  defined  potential  functions.  It  follows  that  the  air  gap 
(which  is  both  current  and  charge  free)  fields  may  be  obtained 
from  the  solution  of 
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The  function  %  is  the  vector  potential  of  the  field  and  can  be 
uniquely  determined  by  satisfying  the  boundary  conditions.  The 
use  of  the  potential  function  A  simplifies  the  exact  solution  of 
the  field  problem  by  cutting  the  algebra  because  all  the  com¬ 
ponents  of  the  field  vectors  are  included  in  it. 

In  cylindrical  co-ordinates  one  may  write  equations  (3.1)  and 
(3.4)  as 


XT  A  -  i  f  V72A  -  2dAe 
V  A  -  lr  ^  V  Ap  vi 


Aj 

T1 


:) 


r  (  72a  +  a  Mr  _ 
v  j  "0  T1  bo  r2 


)1 


.3.5 


+  izv2az  =  Mo 


\7A  -  ■  At  ,  +^AZ. 

VA‘  T  T  +Tie  + 


•  • 


.3.6 


(1)  The  machine  length  has  been  assumed  to  be  considerably  larger 
than  its  radius,  thus  making  the  effects  of  the  currents  flowing 
in  the  end  wires  negligible.  (As  a  matter  of  fact,  some  fringing 
of  the  fields  must  occur  at  the  ends  of  the  machine  and  this  may  be 
neglected  unless  the  length  is  rather  short.) 

(2)  Furthermore,  it  is  now  postulated  that  the  surface  current 
density  is  independent  of  the  z-co-ordinate  of  the  system  (that  is, 
current  is  the  same  in  the  entire  path  in  the  z-direction  and, 
also  there  is  no  retardation  in  the  z-direction.  In  a  real  machine 
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this  is  approximately  equivalent  to  neglecting  intertum  cap¬ 
acitances  and  leakages  between  conductors  of  the  same  member, ) 
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The  above  two  conditions  give 

"b  A 


=  o 


and  taking  A  =  A  = 
r  9 

V.A  =  0  is  satisfied. 


0,  the  divergence  condition,  namely  equation 
The  wave  equation  is  thus  reduced  to 


V  ^A'ToB 
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>t)  =  /<, 


Oeo 


bZA; 


at  9  ,t) 


3  7 

«  s  e  O  o  / 


and  this  has  to  be  solved  for  the  specified  boundary  conditions. 
3.3  Boundary  Conditions 

.  -V  a* 

Since  B  =  VxA  and  A  has  only  a  z-component 

r  -"i  x  bA2  -  /■_  bAz  A 

-  “r  F  «l  A 

The  boundary  conditions  discussed  in  section  (2.4)  may  be 
reduced  to 

(1)  at  the  stator  and  rotor  surfaces 


(r  =  b  and  r  =  a) 


b  A  z 


be 

(2) 

at  the  stator  surface 

b  Az 

br 

T=  b 

and 

at  the  rotor  surface 

b  A  z 

b  r 

r=3L 

is  continuous 
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CHAPTER  4 


SOLUTION  OF  THE  WAVE  EQUATION 
IN  OPERATIONAL  FORM 

4.1  Introduction 

In  general,  the  solution  of  the  three  dimensional  wave  equation 

V 1 A  -  /<■  £  _  n  ...3.1 

7  n  /  13  ^  t£  ^ 

is  needed  for  the  analysis  of  the  machine  especially  under  the 
transient  condition.  However,  for  the  model  chosen,  if  the  current 
density  distribution  is  independent  of  the  z-co-ordinate  of  the 
system  even  under  the  transient  condition,  the  solution  of  equation 


V2Az(r,8>t)=  /*o e0  ^f0r,9,t) 

o  t 


...3.7 


will  be  sufficient  for  the  study  of  machine  performance.  (From  now 
on  A  (r,9,t)  will  be  referred  to  as  A) .  An  equation  similar  to 
(3.7)  is  valid  for  each  of  the  phases  of  the  stator  and  rotor  with 
the  reference  system  being  stationary  with  respect  to  the  winding 
concerned. 

The  wave  equation  may  be  solved  by  various  methods For 

ordinary  steady  state  analysis  of  the  machine,  a  classical  approach 

is  to  seek  cylindrical  wave  functions,  which  satisfy  equation  (3.7) 

subject  to  the  boundary  conditions  of  section  (3.3).  In  such  a 

steady  state  solution,  the  frequency  is  equal  to  that  of  the  driving 

function  (surface  current  density).  In  the  transient  case  the 

initial  conditions  must  be  considered  and  the  wave  equation  may 

42 

be  solved  as  an  eigenvalue  problem  .  Imposing  the  boundary 
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conditions  upon  a  solution  obtained  by  Bernoulli's  separation 
method,  the  eigenvalues  and  the  corresponding  eigenfunctions 
of  the  problem  may  be  determined.  To  satisfy  the  initial 
conditions  it  is  necessary  to  be  able  to  express  an  arbitrary 
function  as  an  infinite  series  of  the  eigenfunctions. 

By  using  the  method  of  the  Laplace  Transform,  an  operational 
solution  of  the  wave  equation  (3.7)  is  shown  in  this  chapter.  The 
general  solution  thus  obtained  can  be  conveniently  applied  for 
certain  analysis  of  the  machine  problems  in  both  the  steady  and 
transient  states  (see  chapters  5,  8  and  9).  Also,  the  general  solu¬ 
tion  may  be  approximated  to  define  frequency  independent  parameter 
for  any  transient  and  steady  state  analysis  of  the  machine  from  a 
conventional  quasi-static  point  of  view  (see  chapter  7). 

4.2  The  Laplacian  Transformation 

The  Laplace  transforms  are  a  pair  of  functions  related  by  two 
fundamental  equations  of  the  Operational  Calculus: 

(a)  the  Laplace  transformation 

co 

res)  =jfct)e.~stat 


and  (b)  by  the  Fourier  -  Mellin  theorem,  the  inverse  formula 

7+> 

f  ct)  -  fres)  £stds  t  >/  0 

y-j® 


The  integrals  are  supposed  to  converge,  with  t  limited  to  real 
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values,  while  s  is,  in  general,  a  complex  number.  The  convergence 
of  the  first  integral  is  assured  if, 
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00  CO 

^  |  f(t)e  s_t | dt  = 

0  ° 

converges.  The  line  integral  in  the  inversion  formula  is  evaluated 
by  transformation  of  the  path  of  integration  to  a  suitable  closed 
contour  and  the  use  of  the  Calculus  of  Residues .  The  inverse  trans¬ 
formation  is  a  solution  of  the  Laplace  transformation  and  vice 
versa,  so  that  the  functions  related  by  one  of  the  above  equations 
are  of  necessity  related  by  the  other  also. 

4.3  The  Operational  Equation  1 
Expand  equation  (3.7)  as 


r, ,  n  |  -Re(s).t,, 
f(t)  e  dt 


V2A  = 


3  2  A 

3  r2 


+  i 


3  A 
r  3  r 
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32A 


2  2 
C  3  t 


. . .4.1 


where 


c  = 


1 


/y  £ 
o^o 


The  solution  of  this  equation  satisfying  two  boundary  and  two  initial 
conditions,  is  sought.  The  initial  conditions  are 
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t->o 


A(r,0,t)  =  A  (r,6)  =  A 
o  5  o 


3  A, 


Lim  g~pr, 0  ,t)  =  A-^(r,0)  =  A-^ 
t-K5 


Take  the  Laplace  transform  of  (4.1)  assuming  that  the  integrals 
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for  s>o,  exist  and  the  operations  of  differentiating  with  respect 
to  the  spatial  co-ordinates  and  taking  the  Laplace  transform  can 
be  interchanged  such  that 
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Therefore,  with  the  above  assumptions  as  to  the  nature  of  the 
vector  potential  A,  one  obtains  from  equation  (4.1)  and  the 
initial  conditions  the  operational  equation  in  the  s -domain  as 
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This  is  the  two  dimensional  form  of  the  Helmholtz  Equation  and  has 
to  be  solved  to  satisfy  the  boundary  conditions. 

4.4  Solution  of  Helmholtz;  Equation 

For  the  complementary  function  write 


p2-  Y1-  cL 


Also,  let 

f(r,  e,s)  =  qCr^Fpe) 


Substitution  of  equation  (4.4)  in  equation  (4.3)  and  a  few 
simple  manipulations  yield 


1  cffl  1  JFi 

fi  drz  rFx  dr 


...4.5 


Let 


From  which, 


Sin  (ing  + 


u. 


6 


(n  being  the  order  of  the  space  harmonic) . 
Substitute  this  in  equation  (4.5)  to  obtain 
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This  is  the  Modified  Bessel  Equation  and  its  solution  is 


F  =  A  I  (-r)  +  B  K  (-r) 

1  n  n  c  n  n  c 

s  s 

where  I  (— r)  and  K  (— r)  denotes  the  Modified  Bessel  functions  of  the 
n  c  n  c 

*  s 

first  and  second  kind  respectively,  both  of  order  n  and  argument  — r. 

Put 


T  t 


-  =  k,  C  =  A  C 
c  n  n  n 


T  t 

and  D  =  B  C 

n  n  n 


Hence,  the  complementary  function  is 


^(r,0,s) 


C  I  (kr)  +  D  K  (kr) 
n  n  n  n 


sin(n6  +  3  ) 
n 


. . .4.8 


where  C  ,D  and  B  are  arbitrary  constants  to  be  determined  from 
n  n  n  J 

the  boundary  conditions . 

For  the  particular  integral  corresponding  to  the  right  hand  side 
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of  equation  (4.2),  Lagrangian  method  of  variation  of  parameters  may 
be  used.  However,  it  is  more  directly  obtained  as  follows: 

Assume  that 

fj  (mx)Y  (mr)  -  Y  (mx)J  (mr)l 
L  n  n  n  n  -l 


(A)  A  (r,0)  =  R 


and 


[~J  (ma)Y  (mb)  -  Y  (ma)J  (mb)| 
L-  n  n  n  n  -1 


sin(n0  +  8  ) 
n 


fj  (mx)Y  (mr)  -  Y  (mx)J  (mr)l 
(B)  A_  (r,  0)  =  R  — £ - 2 - F - F_ - 1 

[j  (ma)Y  (mb)  -  Y  (ma)J  (mb)l 
L  n  n  n  n  J 


sin(n0  +  6  ) 
n 


where,  J  (mr)  and  Y  (mr)  are  the  Bessel  functions  of  the  first  and 
n  n 


second  kind  respectively  of  order  n,RQ,R,m  =  ^  (go  being  angular 
frequency)  are  constants  and 
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x  =  a  or  b  (constants).  In  chapter  8  it  will  be  shown  that  the 
steady  state  solution  of  the  vector  potential  due  to  the  stator 
currents  will  contain  a  factor  similar  to  (A)  with  x  =  a  and,  that 
due  to  the  rotor  currents  will  contain  a  similar  factor  with  x  =  b. 
The  assumption  of  the  particular  form  of  the  initial  conditions 
Ag(r,e)  and  A^(r,e)5  given  by  (A)  and  (B)  respectively,  implies  that 
the  transients  will  be  created  by  a  sudden  change  when  the  machine 
is  operating  at  the  steady  state  with  a  single  frequency  current  in 
the  phase  concerned. 

Let  the  particular  solution  sought  be  of  the  type 

Ap  =  1dn  Sim  (me 

(P  denoting  particular  integral) 

Substitution  into  (4.2)  yields 
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4  4 

2- 

Adding  and  subtracting  m  R  ,  the  above  equation  may  be  reduced  to 
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Since  Y^Cmr)  and  J^Cmr)  are  the  solutions  of  the  Bessel  Equation 
of  order  n, 

fsR  +R,  ~)  fj  (mx)Y  (mr)-Y  (mx)J  (mr)l 
—  _  ^  o  1  7  u  n  n  n  n  -1 

n - 2 — 2 — 77  •  - * - 5 - 1 - B 

c  (k  +m  )  1 J  (ma)Y  (mb)  —  Y  (ma)J  (mb) 

J  n  n  n  n 

...4.11 

The  solution  of  equation  (4.2)  for  the  nth  space  harmonic  is, 
therefore,  given  by 


A  (r,0,s)  =  [c  I  (kr)  +  D  K  (kr)l  sin(n0  +  B  ) 
n  ’  5  Lnn  nn  J  'll 


(sR  +R-, )  fj  (mx)Y  (mr)-Y  (mx)J  (mr)1 
v  o  1*  >-n  n  n  n-1 

— 2 — 2 — 27  *  T1 - ' - 1 - » - T 

c  (k  +m  )  J  (ma.)Y  (mb)-Y  (ma)J  (mb)| 
y  L  n  n  n  n  J 


sin(n0+/6’n) 


...4.12 


The  general  solution  for  A  is  the  summation  of  A  over  n  such  that 


A(r,0,s)  =  )  A  (r,0,s) 

/  i  n 

n=l 


. .  .4.12A 


As  stated  before,  the  initial  conditions  introduced  in  obtain¬ 
ing  equations  (4.12)  and  (4.12A)  correspond  to  a  single  time  frequency 
w.  If  more  than  one  time  function  (conceivably  a  Fourier  series  of 
time  harmonics)  is  involved,  there  would  be  initial  conditions  sim¬ 
ilar  to  (A)  and  (B)  for  each  of  the  time  frequencies  present.  Each 
set  of  the  initial  conditions  would  give  a  corresponding  particular 
integral  and  the  most  general  solution  of  equation  (4.2)  to  include 
all  such  cases  may  be  expressed  as 


A(r,0,s)  = 


A  (r,0,s) 
j  v,n 


...4.12B 
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where,  y  is  the  order  of  the  time  .function  (harmonic  in  the  case  of 
Fourier  Series).  However,  in  this  analysis  only  one  set  of  initial 
conditions  (V  =  1)  will  be  considered. 


4.5  Boundary  Conditions  and  Vector  Potentials  due  to  the 

Phase  Currents 

The  boundary  conditions  in  the  s-domain  are 

(1)  at  the  surface  of  excitation  the  change  in  the  tangential 
component  of  the  magnetic  field  intensity  in  the  s-domain  is 
equal  to  the  linear  current  density  in  the  same  domain. 

(2)  at  the  other  surface  the  tangential  component  of  the 
magnetic  field  intensity  vanishes. 

It  is  convenient  to  use  superposition  and  hence  the  vector  po¬ 
tential  corresponding  to  each  phase  current  will  be  determined. 

(A)  Stator  Phase  (a):  For  the  stator  phase  (a), 


(NOTE:  The  transform  of  the  linear  current  density  (am.p/m)has  been 
replaced  by  that  of  the  actual  current  (amp)  multiplied  by  a  constant. 
See  section  2.5.) 
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s  s 

where,  Rq.,  and  are  constants  and  is  the  angular  frequency  of 
the  stator  currents. 

Using  the  above  relations  with  equation  (4.12)  one  obtains 
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and 
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(B)  Stator  Phase  (b) :  Following  the  same  procedure  as  for  the 
stator  phase  (a) ,  the  vector  potential  for  the  stator  phase  (b) 
is  obtained  such  that, 
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(C)  Rotor  Phase  (a) :  The  boundary  conditions  are 
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where,  R  and  R,  are  constants  and  w  is  the  angular  frequency  of 
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the  rotor  currents . 
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(D)  Rotor  phase  (b) :  Similar  treatment  of  the  rotor  phase  (b) 
gives 
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The  resultant  air  gap  vector  potential  for  all  the  currents 
acting  together  is 
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CHAPTER  5 


ON  THE  INVERSE  TRANSFORM 

5 . 1  Introduction 

The  operational  solution  of  the  wave  equation 

q2.  A 

VzA(T,i.t )=  -frr  ...3.7 


assumed  the  space  co-ordinates  r  and  6  as  independent  of  time 
t  so  that  the  evaluation  of  the  integral  in  taking  the  transform  of 
A  and  operations  of  taking  partial  derivatives  with  respect  to  r  and 
0  are  interchangeable.  In  equations  (4.13)  through  (4.16)  if  the 
current  transforms  are  known,  in  principle  the  inverse  transforms 
may  be  calculated  in  a  routine  manner.  In  general  each  current 
transform  may  be  determined  from  an  equation  similar  to 


E0 


St 
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where  c~is  the  resistivity  of  the  surface  material,  E^  is  the  field 

Aa 

intensity  applied  from  external  source  and  is  the  induced  field. 

Equation  (5.1)  requires  the  stator  vector  potentials  referred 
to  the  rotor  co-ordinate  system  and  vice  versa.  In  determining 
and  its  transform,  motion  of  the  rotor  must  be  taken  into  considera¬ 
tion.  Equations  (4.13)  through  (4.16)  can  be  easily  manipulated  to 
account  for  constant  rotor  speed  at  sinusoidal  steady  state  (see 
chapter  8).  However,  for  the  expressions  to  be  obtained  in  this 
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chapter,  the  rotor  will  be  considered  as  stationary. 
5.2  Electric  Field  Intensities  at  the  Stator  and  the 
Rotor  Surfaces. 

In  the  time  domain 
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Hence,  at  the  stator  surface  (r  =  b) ,  the  nth  harmonic  field 
intensities  are 
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where, 
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Similarly  at  the  rotor  surface  (r  =  a)  the  field  intensities  are 
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5. 3  Line  Integrals 

Note  that  the  resolution  of  the  actual  current  into  Fourier 
space  harmonics  of  linear  current  density  may  be  taken  as  that 
each  of  these  harmonics  is  produced  by  a  fictitious  winding  whose 
conductors  are  sinusoidally  distributed  in  space.  For  every  cycle 
of  the  fundamental  space  component  there  are  n  cycles  of  the  nth 
harmonic  winding.  Let  any  of  these  n  cycles  be  called  the  nth 
harmonic  coil  of  the  nth  harmonic  winding. 

Take  the  Laplace  transform  of  equation  (5.1)  and  integrate  it 
to  obtain 
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where  corresponds  to  the  total  path  of  integration  in  the  nth 


harmonic  coil, 


In  equation  (5.27),  j  R^i  I(S)dfcp  gives  the  total 
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drop  of  the  nth  harmonic  coil,  Eon(s )d(Lp  gives  the  nth  harmonic 
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applied  voltage  and  £  (S)df.p  takes  care  of  the  total  induced  volt- 
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ages.  For  the  last  integral  consider  the  intensities  corresponding 
to  a  conductor  at  net  in  series  with  one  at  not+TT 
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5.1  -  Harmonic  Conductors  in  Series 
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Where  C  is  one  of  the  constants  defined  in  the  previous  section,; N  and 
< are  dummy  variables . 

(A)  Stator  Phase  (a) :  The  conductor  density  is  distributed  as 
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(B)  Stator  Phase. (b):  The  conductor  density  is  distributed 
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(C)  Rotor  Phase  (a) :  The  conductor  density  is  distributed  as 
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N^sin  n(0  -<£>)  and  the  limits  of  integration  are  from  <fi>  to  f+ir/n. 

The  integration  gives 
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(D)  Rotor  Phase  (b) :  Similar  treatment  gives  for  the  phase  (b) 
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Equations  (5.29)  through  (5.32)  are  derived  for  the  nth  harmon¬ 
ic  of  a  2-pole  2-phase  wound  rotor  machine  with  conventional  distri¬ 
buted  windings  on  both  the  members.  Simultaneous  solutions  of  these 
equations  would  give  the  current  transforms,  which  may  be  put  back 
to  equations  (4.13)  through  (4016)  for  the  solutions  of  the  vector 
potentials  in  the  s-domain.  Electromagnetic  angular  momentum  flow  at 
the  starting  of  the  rotor  (rotor  stationary)  may  be  studied  for  any 
transient  by  taking  the  inverse  transforms  of  the  vector  potentials. 
Of  course,  the  transients  must  be  within  the  assumption  made  in 
solving  the  wave  equation. 

For  the  sleeve  rotor  machine  take  <p  -  0  and  write  the  rotor 
equations  point  wise  and  not  in  integrated  form.  This  will  amount 
to  a  change  in  scale  in  the  induced  voltage  terms.  (The  sleeve 
rotor  case  is  made  clear  for  the  fundamental  space  component  in 
chapter  8  and  9 . ) 

5.4  s-Domain  Solution  for  a  Space  Harmonic 

Consider  cp~  0.  This  will  make  the  (a)  phases  independent  of 
the  (b)  phases  and  in  general  for  this  condition  the  solution  will 
be  valid  for  the  sleeve  rotor  machine. 
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Using  the  above  equations  with  (4®  13)  and  (4d5)  one  can  write 
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By  a  similar  procedure  I^Cs)  and  1^  (s)  may  be  obtained  to 


__  g  _ ^ 

determine  A^n(r,0,s)  and  (r,  0,s). 
5 . 5  Discussion 


vr 


The  equations  of  this  chapter  are  suitable  for  the  study  of 
the  machine  performance  characteristics  at  starting*  In  the  last 

_SS  o  o  SS  2?I? 

section  if  V  (s)  is  a  known  function  and  R  and  R  are  known 
an  an  an 

parameters ,  then  in  principle  the  inverse  transforms  may  be  evalua¬ 
ted  immediately.  In  general  these  may  not  be  known  explicitly  and 
interconnections  of  the  harmonic  windings  may  be  necessary  to. obtain 
the  final  solution  in  the  time  domain „  However,  for  a  machine  if 
the  harmonics  are  negligible,  then  the  equations  are  immediately 
applicable  for  the  space  fundamental. 
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CHAPTER  6 


ELECTROMAGNETIC  ANGULAR  MOMENTUM 
FLOW  AND  MACHINE  TORQUE 

6 . 1  Introduction 

As  in  any  electromagnetic  system,  the  law  of  conservation  of  electro¬ 
magnetic  angular  momentum  may  be  conveniently  used  to  determine  the 
production  of  torque  in  a  machine.  With  the  permeability  of  iron 
being  considered  infinitely  larger  than  that  of  air,  the  field  quan¬ 
tities  are  vanishingly  small  inside  the  iron  parts  of  the  machine  and 
mementum  can  only  be  stored  in  the  air  gap.  The  angular  momentum  flow 
balance  is  then  uniquely  given  by  equation  (1.17A)  such  that 


rotor 


where  Vrotor  is  the  volume  of  the  rotor  and  V  is  the  volume  enclosed 
by  the  surface  of  integration.  Of  course,  the  field  quantities  are 
referred  to  the  same  co-ordinate  system. 

If  V  =  V  ,  i.e.  ,  the  surface  of  integration  coincides  with  that 
of  the  rotor. 


Tr  (t)  =  r  x (M.n)dS 


6.1 


rotor 


gives  the  mechanical  torque  on  the  rotor.  The  stator  torque  is  the 
sum  of  this  torque  and  the  rate  of  change  of  angular  mementum  in 
the  air  gap  such  that 
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■p  s 

For  sinusoidal  steady  state  the  time  averages  of  T  and  T 

J  &  em  em 

are  equal  and  is  given  by 
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r  x  (M.n)dS 
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av 


...  6.3 


with  S  being  placed  anywhere  in  the  air  gap. 

6.2  Integration  of  r  x  (M.n)ds. 

From  Maxwell’s  stress  tensor  it  follows  that 


M.n  =  e  (n.M  +  y  (n.8)S  -  n(i  e  E2  +  iy  H2) 
o  o  2o  2  o 

Since  ?!  is  completely  in  the  z-direction,  n.S  =  0  and  the 

contribution  of 
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Furthermore,  it  is  known  that 
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Using  equations  (4.13)  through  (4.16)  or  the  corresponding  time 
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functions  and  performing  the  differentiations  indicated  in  equations 
(6.5)  and  (6.6),  one  obtains 


l3ss  =  i  H  cos  n9s  -  i  H  sin  nes  ...6.7 

an  r  anr  9  ane 

^bn  =  ^rHbnrCOS  n(eS_T7/2)  "  ieHbnesin  n0S  “  T7/2)  ...6.8 

Hrr  =  I  Hrr  cos  n(es-4>)  -  TQHrr  sin  n(BS-<t>)  ...6.9 

an  r  anr  9  ane 


=  i  Hrrcos  n(eS-4>-'n’/2)  -iflHrr  sin  n  (6  S-<t- n/2)  ...6.10 

r  bnr  bn6 


In  particular  cases  the  quantities  may  be  determined  in  s 
or  t  domain  from  the  corresponding  vector  potentials.  Knowing  the 
6  dependence  of  the  magnetic  fields  and  proceeding  with  the  integra¬ 
tion  indicated  in  equation  (6.4)  it  can  be  shown  that  there  would 
be  no  interaction  between  the  magnetic  field  components  of  the  same 
phase.  Furthermore,  there  would  be  no  resultant  contributions  due 
to  the  interaction  between  the  magnetic  fields  of  the  same  member 
and  there  would  be  no  interaction  between  harmonics  of  different 
orders .  Hence 
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Proceeding , 
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Therefore,  the  angular  momentum  stored  in  the  air  gap 

Hu 
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where 


S9  =  CE”  4  Ebn  +  4  KJ  •* 


^Hanr+  H§nr  4  “Sr  4  «"r) 


It  may  be  noted  that  it  would  be  convenient  to  carry  out  the  integra¬ 
tion  in  (6.12)  for  the  particular  cases. 

6.4  Stator  and  Rotor  Torque  Unbalance. 
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The  Stator  Torque 
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The  Torque  Unbalance  is  given  by 
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As  stated  before, 
term  is  zero . 


Rate  of  change  of  stored  angular  momentum 
in  the  air  gap* 

at  sinusoidal  steady  state  the  contribution  of  this 


CHAPTER  7 


SPACE  HARMONICS  UNDER 
QUASI -STATIC  APPROXIMATIONS 

7.1  Introduction 

All  the  equations  so  far  considered  are  in  terms  of  stationary 
co-ordinate  systems  and  the  rotor  has  been  taken  as  yet  stationary. 
When  the  rotor  is  in  motion  the  determination  of  electric  field 
intensities  to  write  the  expressions  for  currents  requires  the  trans¬ 
formation  of  the  vector  potentials  (since  all  other  field  vectors 
may  be  derived  from  the  vector  potentials)  between  the  stator  and 
rotor  co-ordinate  systems,  which  are  in  relative  motion. 

Consider  two  inertial  frames  S  and  S’  in  relative  motion  with  a 
constant  linear  velocity  (relative)  v  along  the  x-axis.  If  $'  is 
the  scalar  potential  and  A'  the  vector  potential  at  a  point  in  the 
reference  system  S’,  it  may  be  shown  that  under 
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Fig.  7.1  -  Inertial  Frames 
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Lorentz  Transformation  the  potentials  transform  as 


<£>  = 


s[T~-  xs'Vc2 


■  (§  -h  19  Ax) 


Ax.  = 


y  /-  -\s>'2/c7- 


(**  +  &*) 


Ay  = 


'y- 


A  <7  -  A  4 


Therefore,  for  a  field  problem  if  =  o,  the  vector 

potential  (with  only  z-component  existing)  remains  invariant  under 
this  transformation  such  that 


Az c x>(y,z,0  =  Az  (x^z.t) 

Let  «c  which  will  make  x=x’  +vt,  t  =  t/  and  the  length  mea¬ 
sures  in  the  two  systems  the  same.  By  applying  a  transformation  of 
the  form 

y/a^ 

B  -  and  T  =  <3_z  £ 

(a^,  a?  and  a^  are  arbitrary  constants) 
the  S  system  may  be  made  to  coincide  with  the  stator  reference  system 
and  the  ST  system  to  the  rotor  reference  system,  which  would  be  in 
uniform  rotation. 

Hence,  it  can  be  seen  that  for  constant  rotor  angular  velocity  w  , 
the  vector  potentials  (z-component  only  exists  in  the  respective 
rest  system)  remain  invariant  under  a  transformation  between  the  stator 
and  the  rotor  co-ordinate  systems.  To  obtain  the  field  vectors  from 
the  partial  derivatives  in  a  system,  stator  or  rotor,  the  vector 
potential  may  be  easily  expressed  in  terms  of  the  co-ordinates  of  that 
system.  In  the  stator-system  if  the  rotor  position  is  defined  as 
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<£  =  i jo  t+<£>  ,  the  co-ordinates  will  transform  as 
w  he 


6r--9s-^J-4, 
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tr=  tS  =  t  •••7-4 

It  may  be  noted  that  this  is  similar  to  the  Galilean  transformation 
of  the  classical  mechanics  and  is  valid  for  the  tangential  velocity 
of  the  field  point  small  in  comparison  with  the  velocity  of  light. 

In  principle  the  operational  solution  of  the  wave  equation 
can  easily  accomodate  the  transformation  of  the  co-ordinates  for 
constant  rotor  speed.  However  if  the  rotor  position  is  given  by 
any  arbitrary  function  of  time,  it  is  convenient  to  use  the  quasi¬ 
static  approximations  of  the  vector  potentials.  The  quasi-static 
approximation,  which  is  in  effect  the  solution  of  Laplace's  equa¬ 
tion  with  the  time  function  written  as  a  factor,  has  the  advantage 
that  frequency  independent  parameters  may  be  derived  (follows  from 
the  frequency  independence  of  the  co-efficients  of  the  time  functions 
involved  in  the  solution  of  the  field  problem)  to  define  the  machine 
equations  (the  volt  ampere  equations  and  the  torque  equation  -  collect¬ 
ively  known  as  the  equations  of  motion) ,  which  automatically  takes 
care  of  the  motion  of  the  rotor. 
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The  solutions  of  these  equations  my  determine  the  machine  char¬ 
acteristics  (transient  and  steady  state).  However,  such  approximation 
my  not  provide  an  adequate  picture  of  the  subtle  effects  taking  place 
in  the  electromagnetic  system  of  the  machine. 

7.2  Quasi-static  Approximations  of  the  Vector  Potentials 

The  quasi-static  approximations  my  be  obtained  by  using  the  low 
argument  approximations  (valid  at  low  frequencies  only)  of  the  Bessel 
functions  involved  in  the  solutions  of  the  vector  potentials.  This 
follows  from  the  fact  that  for  the  frequency  equal  to  zero,  equation 
(4.2)  or  its  counterpart  in  the  time  domin  is  reduced  to  Laplace’s 
equation.  Hence  in  the  limit  as  the  frequency  tends  to  zero  the 
solution  of  either  of  these  equations  will  be  the  solution  of  Laplace’s 
equation  with  the  time  function  or  its  transform  coming  as  a  factor. 
This  is  the  quasi-static  solution  for  the  vector  potential.  Now  for 
z-*0,  the  function1^ 
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In(z)  =  — —  K(z)  =  2n_1(n-l)l  z-n 

2  nj 


n  2n(n-l)  z-n 

Jn(z)  =  —5-  Y  (z)  = - - - 

2n  n!  "  77 

From  the  above  relations  it  my  be  shown  that 
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If  these  relations  are  used  in  equations  (4.13)  through  (4.16),  they 
are  reduced  to 
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The  currents  I_„,  L  ,  It,-  and  I,  will  be  referred  as  harmonic 
an 3  bn 7  an  bn 

currents  and  may  be  taken  as  function  of  time. 


7. 3  Volt-Ampere  Equations 

The  air  gap  electric  field  intensities  are  given  by 
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(Differentiation  in  the  rest  system.) 
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Integrating  these  intensities  as  in  section  (5.3)  one  can  write 
the  volt-ampere  equations  for  the  nth  harmonic  as 
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are  the  nth  harmonic  inductance  parameters  (for  one  coil  of  the 
n  coils).  The  operator  p  stands  for  time  differentiation.  Note 
that  the  inductances  expressed  as  functions  of  the  rotor  position 
which  is  time  dependent.  If  p  denotes  total  differentiation,  the 
effect  of  motion  on  the  electrical  system  is  automatically  included 
in  the  electrical  equations . 

The  equations  (7.13)  are  a  set  of  nonlinear  differential  equa¬ 
tions  ,  the  solution  of  which  is  still  a  great  problem.  For  any 
arbitrary  set  of  conditions  a  total  solution  may  not  be  feasible 
even  if  only  one  harmonic  is  present.  If  the  rotor  position  is  known 
as  an  explicit  function  of  time,  the  equations  are  reduced  to  linear 
differential  equations  with  time  varying  co-efficients  and  may  be 
solved  by  suitable  changes  of  variables  (the  so  called  co-ordinate 
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transformation  of  the  generalized  machine  theory) . 

For  transient  analysis  at  constant  rotor  speed  the  equations 
may  be  referred  to  a  set  of  co-ordinate  systems  which  are  station¬ 
ary  with  respect  to  each  other  and  the  method  of  the  Laplace 
transform  may  be  conveniently  used.  Transients  at  constant  speed 
implies  that  the  moment  of  inertia  of  the  rotor  and  the  attached 
apparatus  is  so  large  that  in  the  limit  there  is  no  change  of 
rotor  speed  during  the  interval  of  time  considered. 

7.4  The  Torque  Equation 

The  magnetic  field  intensities  are 
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Put  the  radial  and  the  tangential  components  of  the  magnetic 
field  intensities  in  equation  (6.11)  to  obtain 
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If  the  rate  of  change  of  electromagnetic  angular  momentum 
storage  is  neglected,  the  electromagnetic  torque  for  the  nth  harmon- 
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The  torque  equation  of  the  machine  is 


e  wtq 
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moment  of  inertia  of  the  rotating  parts. 


rotational  frictional  co-efficient 
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m 


spring  constant 


T  mech  =  Output  mechanical  torque 
7 . 5  Discussion 

Harmonic  analysis  of  a  two-pole,  two-phase  cylindrical  rotor 

machine  under  quasi-static  approximations  is  shown  in  this  chapter 

and  the  method  may  be  easily  extended  to  any  cylindrical  rotor 

machine.  The  results  may  be  suitably  applied  to  any  particular 

case  by  relating  the  harmonic  currents  and  voltages  to  the  actual 

winding  currents  and  terminal  voltages  of  the  physical  device. 

47 

It  may  be  noted  that  White  and  Woodson  have  obtained  these 
results  as  an  extension  of  their  generalized  machine  analysis. 

Their  method  includes  the  solution  of  Laplace's  equation  in  the 
air  gap,  and  the  determination  of  the  inductances  through  energy 
storage.  The  equations  of  motion  are  written  by  using  the  Euler- 
Lagrange  equation  of  the  classical  mechanics.  They  have  also  shown 
how  to  write  the  volt-ampere  equations  through  the  concept  of  rate 
of  charge  of  flux  linkages  and  deduced  the  torque  equation  by  using 
the  principle  of  conservation  of  energy  and  virtual  displacement. 
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As  may  be  summarized  the  method  shown  here  includes  approximation 
of  the  solution  of  the  vector  wave  equation  to  obtain  the  quasi¬ 
static  solutions  for  the  vector  potentials ,  from  which  all  other 
field  quantities  are  derived  by  simple  operations.  The  volt-ampere 
equations  are  written  by  using  the  line  integrals  of  the  electric 
field  intensities.  The  torque  equation  is  obtained  by  using  Maxwell's 
stress  tensor  and  applying  the  principle  of  conservation  of  angular 
momentum  (neglecting  the  effect  of  rate  of  change  of  electromagnetic 
angular  momentum  storage)  in  an  electromagnetic  system.  The  end 
results  are  the  same  as  those  of  White  and  Woodson. 
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CHAPTER  8 


STEADY  STATE  ANGULAR  MOMENTUM 
FLOW  IN  AN  INDUCTION  MOTOR 

8 . 1  Introduction 

The  theory  thus  far  developed  will  now  be  used  for  the  detailed 
analysis  of  a  two-phase,  two-pole  wound  rotor  induction  machine  to 
show  the  electromagnetic  angular  momentum  flow  at  the  steady  state. 
The  cases  of  a  machine  with  a  conducting  sleeve  on  the  rotor  iron 
and  a  similar  machine  without  any  iron  in  the  stator  or  rotor  will 
be  introduced  to  further  enhance  the  concept  of  the  electromagnetic 
angular  momentum  flow.  It  will  be  considered  that  the  stator  excita¬ 
tion  is  balanced  and  the  rotor  is  running  at  constant  angular 
velocity  o^Cc^  =  0  is  a  special  case).  Actually  the  induction  motor 
is  not  a  constant- speed  machine.  However,  in  most  of  the  practical 
machines,  the  mechanical  response  (i.e.  the  rate  of  acceleration  of 
the  rotor)  is  much  slower  than  the  electrical  response;  hence,  a 
constant  speed  constraint  is  legitimate  for  this  kind  of  machines. 

Under  the  assumption  of  balanced  excitation,  let  the  stator 
currents  be  specified  as 

is  (t)  =  IScos(w  t  -  0S)  ...8.1 

a  s  p 

i^  (t)  =  ISsin(w  t  -  0S)  ...8.2 

b  s  p 


A  self  consistent  analysis  will  show  that  for  constant  rotor  speed 
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the  rotor  currents  will  also  be  balanced  and  let  these  currents  be 

l^C-tO  -  Ir Cos  ...8.3 

l^(i)  =  Sin(^{-8p)  ...8.4 

8.2  Inverse  Transforms  of  the  Stator  Vector  Potentials 

Since  only  the  steady  state  field  components  are  required  to 
study  the  electromagnetic  angular  momentum  flow  at  the  steady  state, 
the  evaluation  of  the  inverse  transforms  in  this  section  is  limited 
to  obtain  the  steady  state  stator  vector  potentials.  In  the  inver¬ 
sion  process  these  functions  are  contributed  by  the  poles  of  the 
transforms  of  the  driving  functions  and  may  be  determined  as  follows. 

If  the  initial  conditions  are  zero,  then  for  n  =  1,  equations 
(4.13)  and  (4.14)  are  reduced  to 

=  /io  Nj?  f ' Xi  V,  S)  I*C4)  sin  0s  ...8.5 

and 

A|i(V,bV)  V*0  Ni  -§-XpCt£)I^0  s^(?  -  " /*•)  •  •  .8.6 

—  ss 

where,  X1  (his)  is  given  by 

xllr,s)=  fe' CsrA)  - 1 C O**  3  __>8i7 

1  ’ '  1  Me.) 


The  Laplace  transforms  of  the  currents  in  equations  (8.1)  and  (8.2) 
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Substitution  of  these  into  equations  (8.5)  and  (8.6)  yield 


Aiicr,.fo  =  Csco;l^f'ng|i 

0  •+  COs  j 

...8.10 


ti-S 


=  Isf-)^«S3  Si -  S  MnSpJ 

(b  +  IO3  J 

.  .  .8.11 


Using  the  inversion  theorem 
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In  the  above  inversion  integrals,  the  poles  at  s  =-J^s  will 
determine  the  steady  state  functions.  To  apply  Cauchy's  Residue 
Theorem  for  the  evaluation  of  these  functions,  close  the  contour 
such  that  it  encloses  the  poles  of  the  integrands  at  s  =  ±rj^s 
but  does  not  pass  through  any  other  pole.  Take  the  residues  associa¬ 
ted  with  these  poles. 


Let 
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Using  the  properties  of  Bessel  functions,  it  may  be  shown  that 
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Summing  up  the  above  residues  and  multiplying  by  2tt  j , 
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Similarly  , 


ASbi07  &)  =-/-<o  A^Cbs)  xiSCr’  c°s)  cos  0s  ^  -Bjt) 


. .8.16 


Aiding  equations  (8.15)  and  (8.16)  one  obtains  the  resultant 
steady  state  air  gap  vector  potential  due  to  the  stator  currents 
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8 . 3  Inverse  Transforms  of  the  Rotor  Vector  Potentials 

Setting  the  initial  conditions  equal  to  zero  and  n  =  1,  equa¬ 
tions  (4.15)  and  (4.16)  may  be  simplified  as 


AjiC r.e's)  =  db-ns)  rj©  sin  e* 
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where , 
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Following  exactly  the  same  procedure  as  in  the  previous  section, 
it  may  be  shown  that  for  the  steady  state 
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The  resultant  steady  state  air  gap  vector  potential  due  to  the  rotor 
currents  is 


l£  Oi ' d/c) %  (y  • If  b/e)  -  Yi  O  V/cj 


.  .  .8.24 


(It  may  be  shown  that  for  any  of  the  currents  specified  in 
section  (8.1),  if  there  were  an  nth  space  harmonic,  the  time  domain 
solution  of  the  corresponding  vector  potential  will  be  of  the  form 
shown  in  this  and  the  previous  sections  except  for  a  change  in  N  and 
the  order  of  the  Bessel  functions . ) 

8 . 4  The  Currents 

The  magnitudes  and  phase  angles  of  the  currents  will  now  be 
calculated.  Indeed,  it  will  follow  that  the  assumption  of  the  forms 
of  the  current  expressions  in  section  (8.1)  is  self  consistent. 

(A)  Rotor  Currents 

Using  the  transformations  of  equations  (7.1)  through  (7.4), 


with  respect  to  the  rotor,  equation  (8.17)  may  be  written  as 
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^  f s_s.  Qo^)-Y,'^Wi 

^ °Nl ^jj ^^yf^^l/t)  -Yl (tJ^/c)  X^rY^ 

[sin(9-r-cJr-t-h9pi-k)J 


8.25 


where  C0y  =  -  £0^ 

From  (8.24)  and  (8.25)  the  electric  field  intensities  at  the  rotor 
surface  are 
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Qr^C! (prVs)  -X'  (^-Wc)  Ji  (Ur^  . 

[Ji(-Va/9  YiPr  bAj  — ' ^'(wTa./e)  uppi-  1>/q) 


Qcos  (.  9f-  iJTi -t  9ff )] 


...8.26 


Their  sum  gives  the  total 
such  that, 


[j/( ^r&/c)  YxC^Aj  -Yj  C^r<a/c) 4(C^Q]  ^ 

[  b/c)  -  s'  (^/e)  ^ 

[CoS  "t  ...8.27 

electric  field  intensity  at  the  rotor  surface 


hi  +Eu 


since  the  rotor  is  short  circuited  this  is  the  only  field  intensity 
which  is  to  be  equated  to  the  rotor  resistive  drop  per  unit  length 
along  the  z-direction 

Ej\  -  RjA^CoS 

a  ...8.28 


Integrating  this  equation  as  in  section  (5.3)  one  can  show  that 


s? 


■  .N  ■  ■  : 


/ 


y  j  ■ 


■ 


■■  •  0-v  V  h;-  '•  }  'v'-t 


ST  5  •  J '  ; 


'  :';-v  >  1 


'  ,•  .  . ‘  :  '*•  •■  **  — 

. .  ■  .  ..  ■  ■>«  .1 


••  •  ~ 

' 

. 


V  -  \*r  -xx  ■ 


■ 


/ 


76 


for  the  rotor  phase  (a) 


IrRr  cos  (co  t 
r 


Qr.  ,  rrTr  .  ,  ,  .r, 

0  )  -  co  L  I  simco  t  -  0  ) 
p  r  r  p 


where 


=  a  L  I  sin(co  t  -  0  -  P  ) 

r  r  p  o 


. . .8.29 


R  =  Total  rotor  resistance  per  phase 
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. . .8.30 
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. . . 8 . 31 


The  L  quantities  are  inductances. 
From  (8.29) 
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=  co  L^I5  sin(co  t  -  0S  -  6  ) 
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And  from  this , 
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9?  =  ei+t+'C+nA 


P  1  To  Y 


...8.33 


where , 


=  t^n  — —  =  rotor  impedance  angle 
x  Kx 


...8.34 


Similar  expressions  will  follow  for  Ir  and  6^  if  the  rotor  phase 
(b)  is  considered  giving  self  consistent  results. 

(B)  Stator  Currents 

With  respect  to  the  stator,  the  rotor  vector  potential  is 
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[sin  Ce-^t  +6p 


...8.35 


From  (8.17)  and  (8.35),  the  electric  field  intensities  at  the  stator 
surface  r  =  b,  are 
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Let  the  applied  voltages  be 

V^Cfc)  =  Vscos  ^st  ...8.38 

\S^60  =  V  s5in  cbst  ...8.39 

Integrating  the  field  intensities  for  the  stator  phase  (a) ,  it  may 
be  written  that 


R^cbsc^t-eb  =  Efjdfp  +  Vscosust 
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R?I scos(Uit-^)-  Is(irn(tos{-  el_) 

=  ^  X*"  ^1^4  (^t _  6p  ^  Cos L0$t  ...8.40 


where . 


Ru 


=  Total  stator  resistance  per  phase 
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Again  the  L  quantities  are  inductances. 
Equation  (8.40)  may  be  solved  as  follows: 
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0^+j^si^)lsCoS  C^t-6p) 


SY'  Y  .  v.  s 

=  LOsL  I  Sin  CHst-Pp  +  £)+V  CoscOst 


Let  - 


where 


gVj4lis=  |^|eJ^ 


lzsl  =  and  ^  -  u rC1  ^sL 

Rs 


SS 


Then  with  exponential  notations 


|ZS|  e^Isej9pS 


-  OJ, 


L^Trti(<Po-^rr/^  ^ 


From  which 


So  that 


l£e1#P=  -i 


IZslL 


0 


-J'fe 


Calculation  of  Is 


Since  |e =  l,  and 


Hence 


Vs-  cos  Lsri  V 

Vs  -  LJsLS*IrCos  C%  -  &p  +  n/g)  rr/z) 


..  .8.41 


Calculation  of  6s 
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It  follows  that 


=  t  s.n 


.1  [Vs  sin^s  -  coaL«*I*Sin  C  rs-4>0  +  &p  -  tt/Q] 

[V  sCoS  ts  ~  CoS  C% -p0  +  0^  -  rr/2)] 


.. .8.42 


Similarly 


as  assumed. 


l^CtD  =li  si^(_h)5f-e*) 


8 . 5  R  otor  Current  Density  of  the  Sleeve  R  otor  Machine 
The  current  density  may  be  calculated  as  follows: 

In  equation  (8.25)  let 

I  sn/=  Xi  a1*1  Bp~t%  =  ^ 

So  that  with  exponential  notation  (the  imaginary  part  gives  the 
solution) 
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By  analogy  with  section  (8.2),  it  follows  that  this  vector 
potential  may  be  produced  by  two  equivalent  current  sheets  placed 
in  space  and  time  quadrature,  the  space  being  coincident  with  the 
stator  but  moving  at  an  angular  velocity  with  respect  to  the 
stator  co-ordinate  system.  In  this  space  the  resultant  current 
density  distribution  producing  A^C^jQ  ,t)  is 

-L 


TS  s  JCOr-^yrt) 


. . .8.44 


Since  the  rotor  is  stationary  with  respect  to  co-ordinate  system 
(-r,9r,t),  the  rotor  current  density  will  have  the  frequency  w  ,  but 
differ  in  magnitude  and  phase  from  the  current  density  in  equation 
(8.44).  Let  the  resultant  rotor  current  density  be 


IiO>t)=  11  QJ 


'C  JC  ~  +  Bp  ) 


...8.45 


The  vector  potential  due  to  this  current  density  is 


b/c)  Y jYY r ~Y-i  3(8-2  -6Jp t  +  9p) 

*A/fe)Yi  C^r  -  %.  (P^/e) 
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To  evaluate  l"T  and  0T,  equate  the  electric  field  intensities 
1  P 

*  on .the  rotor  surface  caused  by  the  stator  and  the  rotor  currents 
to  the  ohmic  voltage  drop  per  unit  length  along  the  rotor  such  that 
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where  R  =  rotor  surface  resistivity. 

r 

Equating  the  real  and  imaginary  parts  of  (8.47),  1^  and  9^ 
may  easily  be  calculated.  Also  referring  the  vector  potentials  to 
the  stationary  stator  co-ordinate  system,  the  current  density  1^ 
and  its  phase  angle  may  be  calculated  in  terms  of  terminal  volt¬ 
ages.  Since  the  form  assumed  for  the  rotor  current  density  in 
(8.45)  leads  to  self  consistent  result,  the  assumption  is  justified. 
8.6  Angular  Momentum  Flow 

Using  equations  (8.25)  and  (8.24)  one  can  write  with  respect 
to  rotor  co-ordinate  system 
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The  time  average  of 


rx(M.n)dS  will  be  evaluated  with  S  being  a 
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However , 


and 
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(The  above  s implication  has  used  the  Wronskian  relation 
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Hence 
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•pc 

where  L  is  given  by  (8.31). 

Equation  (8.50)  gives  the  average  rate  at  which  electromag¬ 
netic  angular  momentum  is  crossing  the  air  gap  of  the  machine  and  is 
absorbed  on  the  rotor  surface.  Since  the  time  average  of 
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is  zero,  the  average  mechanical  torque  (known  as  the  electromag¬ 
netic  torque)  on  the  rotor  is 


Expressing  the  field  components  in  terms  of  the  stator  co-ordinate 
system,  it  may  be  readily  shown  that  the  magnitude  of  the  reaction 
torque  on  the  stator  is 


S  av 

(the  surface  S  being  stationary  with  respect  to  the  stator) 


But 


Tr  -  >pS 
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Hence 


showing  that  the  rate  at  which  the  electromagnetic  angular  momentum 
is  crossing  the  air  gap  is  invariant  under  a  transformation  of  the 
co-ordinates . 

8 . 7  Power  Flow 

The  average  power  flow  across  the  air  gap  is  given  by 


Pa.v  =-  [(EXH).rias  Sr<30 
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with  respect  to  rotor  co-ordinate  system 
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where,, 


=/<oCJS*Ir 
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Some  of  the  terms  in  equation  (8.54)  produce  double  frequency 
components,  the  time  average  of  which  is  zero.  The  net  contribution 
of  the  integral  in  equation  (8.53)  comes  from  E'f^tfyF  and  pcrprs 
with  respect  to  the  rotor  co-ordinate  system,  the  value  of  the  in¬ 
tegral  is 
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It  may  be  readily  shown  that  with  respect  to  the  stator  system 


^v=HLsrIslrsmC6?-B|-t 


...8.58 

Equation  (8.57)  and  (8.58)  reveal  that  the  power  crossing  the 
air  gap  depends  on  the  co-ordinate  system  such  that 
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This  difference  arises  due  to  the  fact  that  with  respect  to  the 
rotor  system,  the  machine  is  not  doing  any  mechanical  work;  -  1  cor 
is  lost  in  the  ohmic  resistance.  Whereas,  with  respect  to  the  stator 
system  the  machine  is  doing  mechanical  work  at  a  rate  Toj^ in  addition 
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to  supplying  the  ohmic  loss.  It  may  be  noted  that  P®  is  the  true 
power  crossing  the  air  gap.  It  comes  from  the  stator  which  is  supp¬ 
lied  from  the  connected  power  system. 

8 . 8  P-T  Relation 

From  equation  (8.51),  (8.52),  (8.57)  and  (8.58)  it  follows  that 
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In  general  it  may  be  written  that 
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Where,  co  is  the  angular  velocity  of  the  field  components  relative 
to  a  co-ordinate  system  with  respect  to  which  the  surface  of 
integration  to  obtain  P  is  stationary,  co  may  not  be  equal  to  the 
angular  frequency  of  the  time  functions  Involved.  For  example,  if 
the  surface  current  distributions  are  sinusoidal  functions  of  n0 , 
instead  of  0 ,  then 
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60, 
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71  m  ^  71 
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Where,  is  the  angular  velocity  of  the  field  components  with  res- 
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co. 


pect  to  the  stator  and  is  the  angular  velocity  with  respect  to 

n 

rotor.  Evaluating  the  integrals  for  P  and  T  it  may  be  shown  that 
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Finally ,  it  may  be  noted  that  the  relation  involving  P  £  and  co  is 
obtained  strictly  from  the  concept  of  electromagnetic  angular  mom¬ 
entum  flow  in  the  air  gap  of  the  machine. 
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8 . 9  Electromagnetic  Angular  Momentum  Velocity  Across  the  Air 

Gap  of  the  Induction  Machine 

A  quantity  of  some  interest  is  the  velocity  of  propagation  of 
electromagnetic  angular  momentum  across  the  air  gap  of  the  induction 
machine.  This  also  requires  a  careful  definition.  By  analogy  with 
fluid  flow,  the  electromagnetic  angular  momentum  velocity  may  be 
defined  as  the  ratio  of  electromagnetic  angular  momentum  flux  per 
unit  area  per  unit  time,  divided  by  the  electromagnetic  angular 
momentum  density.  In  the  case  of  the  induction  motor  with  a  short  air 
gap  length,  time  average  values  of  these  quantities  are  used  and  the 
electromagnetic  angular  momentum  velocity  is  taken  as  the  ratio  of 
the  time  average  electromagnetic  angular  momentum  flux  across  the 
air  gap  to  the  mean  stored  electromagnetic  angular  momentum  per 
unit  air  gap  length,  that  is, 
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Where  the  integration  in  the  numerator  is  over  a  cylindrical  surface 
in  the  air  gap,  g  is  the  air  gap  length  and  the  integration  in  the 
denominator  is  over  the  entire  volume  of  the  air  gap. 

From  equation  (8.52),  at  the  steady  state 
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with 
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For  low  frequencies  and  small  air  gap  length 
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From  equation  (6.12)  the  time  average  of  electromagnetic 
angular  momentum  storage  is 
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With  respect  to  the  stator  co-ordinate  system 
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Expressing  the  field  components  in  terms  of  the  stator  system, 
S0  may  be  obtained  from  equation  (8.64).  Substitution  of  SQ  in 
equation  (8.63)  and  the  evalution  of  the  indicated  integral  gives 
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and 
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Substituting  equations  (8.62)  and  (8.68)  in  equation  (8.61) 
one  obtains 
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where  X  ,  j  ,5  and  5  are  given  by  equations  (8.64),  (8.69),  (8.70) 
and  (8.71)  respectively . 

For  low  frequencies  and  short  air  gap 
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where  ^  and  d  are  given  by  equations  (8.75),  (8.76),  (8.77) 
and  (8.78)  respectively. 

8.10  Energy  Velocity 


As  in  the  case  of  electromagnetic  angular  momentum,  one  might 


enquire  as  to  the  velocity  of  propagation  of  electromagnetic  energy 
across  the  air  gap  of  the  induction  machine.  For  a  short  air  gap 
length  a  definition  of  the  energy  velocity  u>£ ,  similar  to  that  of  the 
angular  momentum  velocity^  ,  gives 
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Where  the  integral  in  the  numerator  is  the  time  average  electromag¬ 
netic  power  flux  across  the  air  gap  and  the  integral  in  the  denomina¬ 
tor  is  the  time  average  electromagnetic  energy  stored  in  the  entire 
volume  of  the  air  gap. 

The  steady  state  power  flux  is  obtained  from  equation  (8.58) 
such  that 
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If  Wg  is  the  average  electric  energy  stored,  it  may  be  shown  that 
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which  is  equation  (8.68)  multiplied  by  o  .  For  low  frequencies  and 
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where  the  quanties  ?  ,  t  £  and  d  are  given  by  equations  (8.75), 

o  1  -d. 5  ->1 

(8.76),  (8.77)  and  (8.78)  respectively. 

If  WMis  the  average  magnetic  energy  stored,  it  may  be  readily 
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where 
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For  low  frequencies  and  short  air  gap  length 
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Where,  P®v,  W£  and  are  defined  by  equations  (8.81),  (8.82)  and 
(8.84)  respectively. 

For  low  frequencies  W  may  be  neglected  in  comparison  with  W 
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and  =  a(b2-a2)  ...8.100 

8.11  Discussion  and  Numerical  Results 

(A)  Discussion  on  the  Exact  Solution:-  A  detailed  field  theory 
analysis  of  an  induction  motor  has  been  developed  to  study  the 
steady  state  production  of  torque  in  the  rotor  in  terms  of  transfer 
of  electromagnetic  angular  momentum  across  the  air  gap.  For 
specified  surface  current  distributions,  the  air  gap  field  quantities 
are  obtained  from  the  exact  solution  of  Maxwell’s  equations.  It  may 
be  noted  that  these  fields  have  the  form  of  electromagnetic  waves. 

An  examination  of  the  nature  of  these  waves  may  be  helpful  to 
clarify  the  physical  implications  involved  in  the  analysis.  Since 
the  general  description  of  the  family  of  waves  remains  the  same 
whether  they  are  referred  to  the  stator  or  rotor  co-ordinate  sys¬ 
tem,  the  discussion  will  be  limited  to  the  latter  case. 

Consider  equations  (8.27)  and  (8.48).  These  are  the  electro¬ 
magnetic  fields  produced  by  the  stator  current  acting  alone.  It 
is  interesting  to  note  that  a  point  (r,0,z)  the  components  of 
the  magnetic  field  intensity  form  an  elliptically  polarized  vector. 

In  a  plane  z  =  constant,  the  electromagnetic  fields  result 
in  a  circularly  polarized  wave.  In  the  three  dimensional  space  of 
the  air  gap  the  field  pattern  is  a  circularly  polarized  standing 
surface  wave.  Being  a  standing  wave  it  gives  storage  of  energy 
without  propagation  in  any  direction.  It  may  be  seen  that  due  to  the 
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interaction  of  the  z-directed  E-field  and  the  radial  component  of  the 
H-field  of  the  circularly  polarized  standing  wave  there  can  be  a  net 
electromagnetic  angular  momentum  storage  in  the  air  gap. 

Gbnsider  equations  (8.26)  and  (8.49),  which  give  the  field  quan¬ 
tities  produced  by  the  rotor  current  acting  alone.  The  field  pattern 
is  of  an  identical  form,  namely  a  circularly  polarized  standing 
surface  wave,  as  that  produced  by  the  stator  current  except  for  a 
phase  displacement  and  a  difference  in  magnitude.  Acting  alone  this 
wave  gives  a  storage  of  electromagnetic  energy  and  angular  momentum 
without  any  propagation. 

Consider  now  the  combination  of  the  stator-and  rotor-induced 
circularly  polarized  waves.  The  interaction  between  the  E-field  of 
each  of  these  waves  and  the  tangential  component  of  the  H-field  of 
the  other  results  in  a  net  flow  of  energy  radially  inward.  Similarly, 
the  radial  component  of  the  H-field  of  each  wave  reacts  with  the 
tangential  component,  of  the  H-field  of  the  other  to  give  a  net  el¬ 
ectromagnetic  angular  momentum  flow.  The  electromagnetic  energy  and 
angular  momentum  are  absorbed  on  the  rotor  surface  at  continuous 
rates  so  that  the  flow  of  these  quantities  is  maintained 

The  concept  of  electromagnetic  angular  momentum  flow  can  be 
further  enhanced  by  considering  a  machine  without  any  iron  but  other¬ 
wise  the  same  as  that  discussed  so  far.  It  may  be  shown  that  the 
field  pattern  produced  by  the  stator  or  rotor  current  is  a  circularly 
polarized  standing  wave  inside  it  (stator  or  rotor)  and  a  radiation 
field  outside  with  waves  spiralling  outwards.  Thus,  inside  the 
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rotor  there  is  a  circularly  polarized  standing  wave  made  up  of 
contributions  from  the  stator  and  rotor  currents .  This  wave  gives 
no  propagation  of  electromagnetic  energy  and  angular  momentum  in 
any  direction. 

In  the  air  gap  between  the  stator  and  rotor  there  is  a  combina¬ 
tion  of  the  stator- induced  standing  wave  and  the  rotor- induced  out¬ 
ward  spiralling  propagating  wave.  The  combined  field  cause  a  net 
flow  of  energy  and  electromagnetic  angular  momentum  inward.  Outside 
the  stator  there  is  a  resulting  outgoing  spiralling  wave  made  up  of 
the  sum  of  the  component  waves  produced  by  the  stator  and  rotor,  An 
interesting  fact  is  that 
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does  not  vanish  when  taken  over  a  surface  just  outside  the  stator 
and  as  a  result  the  mechanical  torque  on  the  rotor  is  not  exactly 
equal  to  the  mechanical  torque  on  the  stator.  The  reason  for  this  is 
that  the  stator  is  radiating  energy  and  electromagnetic  angular 
momentum  outwards  into  space.  In  the  actual  machine  the  presence 
of  an  iron  circuit  outside  the  stator  current  sheet  prevents  (when 
y^)  the  outward  flow  of  electromagnetic  angular  momentum,  thus 
making  the  stator  and  rotor  torque  equal  in  magnitude.  The  main 
effect  of  iron  in  the  machine  is  an  increase  of  the  air  gap  magnetic 
flux  by  a  large  factor,  thereby  increasing  the  torque  and  power 
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(B)  Discussion  on  Quasi-static  Approximations:-  As  noted  in 
Chapter  7,  the  quasi-static  field  quantities  may  be  directly  obtained 
by  using  the  small  argument  approximations  of  the  Bessel  functions. 

At  power  frequencies ,  even  taking  the  mean  radius  of  the  machine 
as  great  as  1  meter,  the  argument  is  of  the  order  of  10  .  From 
the  series  expansions  of  the  Bessel  functions  it  may  be  shown  that 
the  approximations 


J  (z) 
Y1(z) 


give  the  values  of  these  functions  and  their  derivatives  correct  to 
one  part  in  at  least  10^.  Hence,  it  may  be  stated  that  in  mag¬ 
nitude  the  quasi- static  field  values  are  correct  in  about  the 
same  proportion. 

It  may  also  be  noted  that  the  momentum  concept  may  be  used 
to  evaluate  the  torque  even  if  the  electromagnetic  field  quantities 
are  calculated  on  the  basis  of  quasi-static  approximation.  How¬ 
ever,  in  the  ironless  machine  with  such  an  approximation 

rx(M.n)d,S 

Vi 

vanishes  outside  the  stator  showing  no  radiation  of  energy  and 
electromagnetic  angular  momentum  outwards  from  the  machine. 

This  truly  remarkable  agreement  between  the  exact  and  approx¬ 
imate  solutions  only  refers  to  the  steady  state  operation  of  the 
machine.  It  will  be  shown  in  the  next  chapter  that  if  a  sudden 
change  is  made,  for  example  in  the  value  of  the  rotor  resistance, 
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the  quasi-static  approximation  is  inadequate  to  discribe  the  forces 
which  operate  during  the  transition  from  one  steady  state  to 
another.  A  change  in  the  amount  of  electromagnetic  angular  mom¬ 
entum  stored  in  the  air  gap  of  the  machine  must  take  place  and 
there  will  be  torque  impulses  on  the  stator  and  rotor  windings 

to  supply  this  change. 

(B)  Numerical  Results 

The  steady  state  machine  characteristics  such  as  torque- 
slip  relationships  with  rotor  impedance  control  may  be  deLermined 
by  simple  manipulation  of  the  relevent  equations  obtained  in 
this  chapter.  The  calculated  results  should  agree  with  experi¬ 
mental  ones.  Numerical  calculations  are  made  for  a  sleeve  rotor 
machine  with  two  rotor  resistance  values,  to  study  energy  and 
angular  momentum  velocities  as  functions  of  air  gap  length.  The 
rotor  is  taken  as  stationary  and  a  change  in  rotor  surface  res¬ 
istivity  implies  a  change  in  sleeve  thickness . 

Specifications 
a  =  lm  =  Rotor  Radius 

Rr  =  1.72C10"5)  or  1.72C10"4)  ohm/m2 

=  Rotor  Surface  Resistivity 

vll  =  velocity  for  Rp  =  1.72(10  )  ohm/m2 

v^2  =  velocity  for  R^  =  1.72(10  4)  ohm/m2 

The  results  are  shown  as  graphs  in  Fig. (8.1) 

It  may  be  noticed  that  across  the  air  gap  of  an  induction 
machine  the  average  electromagnetic  momentum  velocity  is  higher 
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Fig.  8.1  -  Energy  and  Electromagnetic  Angular  Momentum 
Velocities  Versus  Air  Gap  Length 
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than  the  average  energy  velocity.  The  main  reasons  are 

(1)  The  equations  (8.79)  or  (8.79A)  and  (8.97)  or  (8.97A)  giving 
these  velocities  do  not  account  for  all  electromagnetic  field 
quantities . 

(2)  Because  of  the  term  1/c  the  electromagnetic  angular  momentum 
density  is  a  much  smaller  quantity  than  the  energy  density. 

For  the  same  air  gap  length  both  v^and  v^  have  shown  increased 

values  with  increasing  rotor  surface  resistivity.  The  physical  impli¬ 
cations  may  be  that  for  higher  resistance  values  energy  and  electro¬ 
magnetic  angular  momentum  are  absorbed  at  higher  rates,  so  they  travel 
with  faster  velocities. 
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CHAPTER  9 


ANGULAR  MOMENTUM  FLOW 
AT  TRANSIENT  STATE 


9.1  Introduction 

When  a  sudden  change  is  made  in  a  system,  the  solution  of  the 
system  differential  equations  must  satisfy  the  initial  conditions. 
The  solution  by  the  method  of  Laplace  Transform  takes  care  of  that 
by  finding  out  the  free  response  consisting  of  all  the  terms  (known 
as  transient  terms)  which  are  contributed  by  the  natural  modes,  i.e. 
the  roots  of  the  characteristic  equation  together  with  the  forced 
response  consisting  of  the  terms  contributed  by  the  poles  of  the 
driving  functions.  (The  latter  terms  are  known  as  the  steady  state 
solutions  when  the  driving  functions  are  periodic.) 

The  nature  of  electromagnetic  angular  momentum  flow  and  the 
method  of  its  study  remains  basically  the  same  no  matter  how  the 
transient  is  instigated.  The  problem  will  be  illustrated  by  using 
a  simple  transient.  The  transient  will  be  created  by  assuming  that 
a  two-pole  (space  fundamental  only  present),  two-phase  machine  is 
operating  at  steady  state  with  rotor  open  circuited  and  constrained 
to  remain  stationary.  At  t=0  a  sudden  change  is  made  in  the  value 
of  the  rotor  resistance  (open  circuit  to  closed  circuit).  The 
stator  is  supplied  from  a  "Constant  Current  Source".  If  the 
stator  is  supplied  from  a  "Constant  Voltage  Source"  the  solution 
of  the  problem  will  involve  only  more  algebraic  manipulations. 


1 

^  1  •' 


■ 


' 


v  4 


5  > 


j 


■ 


-  -  '■  ;  ■  .  -  '  •  ■'•  ■  -■  1  •• 'V-:  -  • 

' 

'  ■  '  ,  ■■  ■  ■  '■  ■'  ■  ’  1  -■  1  ■  - 

r 

■ 

, 


■ 


"  •  ,  -*r  ^  ■ 


rci.oiy 


i 

. 

-  .  ..  ■‘1  / 

*  ..  :  .y  .  ■  ■  sk  - 

..  ■■  ;<  ■ '  l  :  '  '  ■■  ’  ■'•  *  •  ; 

.  '  '■  ■  "  >  '■  v  ■ 

- 


102 


9 . 2  Initial  Conditions  and  Resultant  Vector  Potentials 

Let  the  stator  surface  current  densities  (amp/ length)  corres¬ 
ponding  to  the  currents  (amp)  supplied  from  the  "Constant  Current 
Source”  be  specified  as 

1  j  a_(t)  =  if  Cos  «5t 


Ifsin^t  ...9.2 

The  ideal  source  is  taken  to  have  an  infinite  impedance,  hence,  these 
current  densities  will  be  maintained  during  the  transient  conditions 
and  the  corresponding  vector  potentials  may  be  obtained  from  equa¬ 
tions  (8.15)  and  (8.16)  by  substituting  I®  for  N^IS  and  putting 

0s  =  0  such  that 
P 

cr^s3si7i9coswst  •*-9-3 


Ai  i  ^ ^  =-  h  zi  xf  &  °s)  cos  a  Sln  ^ 

where 

yss  Y -YiC^afe>3iC ^/cj] 

^  3  '  [  Ij  OVtyc)  “  'f1C'Wc) 


A2f(r , 0 ,t )  and  A^(r,0,t)  are  known  functions  in  the  time  domain. 

Initially  the  vector  potential  and  its  derivatives  arising  from 
the  rotor  current  are  zero.  Furthermore,  as  implied  in  section  (8.5), 
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the  electric  field  intensities  at  the  rotor  surface  due  to  stator  cur¬ 
rents  vary  as  sine  and  sin(  6-77/2)  and  these  are  the  only  external 
driving  functions  for  the  rotor  current.  Hence,  for  the  same  co-ordi¬ 
nate  system  a  component  of  rotor  current  will  vary  as  sin 6  and  the 
other  as  sin(  8-tt/2)  .  Denoting  the  current  densities  corresponding 
to  the  above  two  components  of  rotor  current  as  1-^  and  1^  and 
setting  the  initial  conditions  zero,  the  s-domain  equations  of  the 
rotor  vector  potentials  may  be  written  from  equations  (4.15)  and  (4.16) 
as 


Arn(r,6,s)  =  y  1^  (s)  —  (r,s)  sine  ...9.5 

al  ’  ’  o  la  si  5 

A^(r,6,s)  =  y  I^(s)  ^  X^r(r,s)sin( 6  -77/2)  ...9.6 

where 

[K  , (sb/c)In (sr/c)-I  ,  (sb/c)K-.  (sr/c)  1 
r^rr,  >.  L  1  1  1  1  4 

X-,  (r,s)  =  —  j - ? - , - , - 

[_K  ^(sa/c)!  ^(sb/c)-!  1(sa/c)K  ^(sb/c^  ...8.20 

As  in  section  (8.5)  I-^(s)  and  I-^(s)  may  be  evaluated  by  equating 
the  resistive  drop  along  the  rotor  length  for  each  of  this  current  to 
the  electric  intensity  Ez(s)  which  is  the  sum  of  the  intensities 
produced  by  the  stator  and  rotor  currents  of  similar  phases  such  that 


R  1^  (s)sin6  =  ESn(s) 
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where 


R  =  El  ©  +  bhcs) 


lY  It 


hi' 


'hi 


...9.8 


R  j  - 


Rotor  Surface  Resistivity, 


E^iCS)  - 


o^Llj 

(SV  ) 


SlYl  9 


...9.9 


T  f c\ ' — 


SCOcLlf  .  , 

-  — ^ — 7T-  ^lnCB'0/2^ 


sVcoJ" 


. .  . 9 . 10 


with 


e2  2_ 


/A<0  O*-  773 


1 


LPi  i^^/cYi^/e)  -Y^/c)^yd] 


.. .  .9.11 


and 


f 


,  .  wirr. 


aloo  =  -/^oc  Iia_(s;  K-  (s;  Sin  9 


. . .9.12 


with 


Hence , 


E11  =  _/M°C  ^ihCS)  *1  ^  S1T1  ^  ...9.13 


LKCsVc)  TgC^J-IiC^V^KxC&^fcjl 
1  ^sa-/e)  Ia  (s  1/q)  -  I^Cs^/c)  Kjsb/c)] 


Ji*Cs)  = 


co; 
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[Rr  +A>c  XxTrCs)J  (  sV  coy 


...9.15 
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T>- 


S^Llf 


(S2+  col )  [  Rr  +  fi0c  %l  r(ji)} 


...9.16 


Putting  these  current  densities  in  the  vector  potential  equations 
(9.5)  and  (9.6)  one  obtains 


V  (X  a  S)=  A  £  L  X1  XIrC^S) 

41  ’  )-t°  S  CsU^j[Rr  +Ac  x«(^] 


.  .  .9.17 


^£iCr>e.s)=Ac 


c  —  Yt  v 

cogLIi  Xx  (t^sJcosq 
CsV^;[Rr+-  x"[rC5i] 


• . . 9 . 18 


9 . 3  Time  Solution  for  the  Vector  Potentials 


The  time  functions  may  be  obtained  from  the  Inversion  Theorem, 


r-y 


CO 


Acu(:r’9>t;)= 


Co  sin  9 

2.77J 


st 


r-j^ 


.  scsVccg-;Q?r+/Aocx«cs51 


...9.19 


where 


do  Cos  9 


y+jco 


xfr(r,ye.stds 


. . .9.20 


2nJ  AJ sn^tRj.+^cxyo] 


d  =  A>cHlIi 


. .  .9.21 


The  integrands  are  single  valued  functions  of  s .  Near  the  origin 
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the  functions 


(r,s)  and 


(s) 


behave  as 


Lim  X  rr(r,S) 
S+o  1 


sa  [r/b  +  b/r] 
(b/a  -  a/bl 


Lim  X^Cs) 
S"K) 


sa  [a/b  +  b/a] 
c  [b/a  -  a/b] 


Hence  the  integrands  do  not  have  any  pole  at  the  origin.  The 
poles  are  at 


and  at 


i  03v 


(corresponds  to  the  driving  function) 


(corresponds  to  the  transient  terms) 


The  latter  poles  are  the  roots  of  R  +  y  c  X-,  (s)=o 

^  r  o  1 

and  these  are  calculated  in  Appendix  ( 3 ) . 

The  first  order  approximation  shown  in  equations  (A3. 38)  and  (A3. 34) 

gives 

T  =  „  =  [a/b  +  b/a] 
r  yo  [b/a  -  a/bj  q 
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2(b-a) 


In 
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. .  .9.23 


nc  Cl  +  2v] 
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v=o,l  ,2 , . . . 


...9.24 
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Fig.  9.1  -  Contour  of  Integration 
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Consider  the  integrals  taken  over  the  closed  contour  of  Fig.  9.1. 
The  circle  of  radius  R  =  does  not  pass  through  any  pole  of  the 
integrands .  In  appendix  ( 4 ) ,  it  is  proved  that  as  v  -*■  °° ,  the  inte¬ 
grals  over  the  arc  FGE  tends  to  zero.  Thus  in  the  limit  as  v  -> 
the  line  integrals  may  be  replaced  by  the  closed  integrals  over  the 
contour  EFGE  such  that 


(r,6 ,t) 


QsinB 


X^r(r,s)esi:  ds 

(s2+co  2)  f~R  +y  cX^  (s)] 
s  L  r  o  1  J 


.9.25 


A^-j^Cr,  6  ,t)  = 


£  cos  e  r 

°  O*- 


2H] 


X^Tr,s)es't:  ds 


(s  +o)  )  Tr  fi  cX^T  (s)l  ...9.26 

s  '  L  r  o  l  J 


Cauchy’s  Residue  Theorem  may  now  be  used  to  evaluate  the  integrals 
of  the  above  equations. 

J? 

Evaluation  of  A^(r,6,t) 

(A)  The  residue 

at  s  =  ico 
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J  s 


.„rr,  ^  -id)  t 
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2n3 
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Q)sin0  ws 
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where 


X,  (r,to  )  is  given  by  equation  (8.23)  with  oo  replaced  by  to  ;  X^ 1  (to  ) 

-L  S  10  S  _L  S 


,rr. 


is  given  by  the  same  equation  with  to  replaced  by  to  and  r  replaced 

T  S 


by  a. 


(B)  The  residues  at  s  =  -a  +i3 

v  —  v 


The  first  order  approximation  will  be  used.  From  appendices  (3)  and 
(4)  it  follows  that  for  large  values  of  s 


Xf  (r,s) 


cosh  [(b-r)s/  ] 
c 


rr 


[R  cXn 
r  Ko  1 


: _ =  &  .  _  i  _  . . 

(  v-i  /r  R  sinhL(b-a)s/c  +  y  coshL(b-a)s/cJ 
(s)J  r  o 


Also  for  large  values  of  s 


,  2  ,  2x  2 

s(s  +  to  )  =  s 
s 


Hence  the  integrand 


77-rr,  x  st 
to  X,  (r,s)e 
s  1  5 


to  cosh  [(b-a)s/c]e 
s 


st 


s(s2 3 *  +  to2)  [R  +  vi  cXf^s)] 
s  r  o  1 


r  3 


s  [R  sinh[(b-a)s/c]  +  y  c  cosh[(b-a)s/c. 
r  o 


The  residue  will  be  evaluated  by  writing  the  integrand  as  a  ratio  p(s)/q  (sq) 
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s~t 

oo  cosh[b-r)s/c]e 
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and  the  series  terms  varying  as  l/#.3  in  equation  (9.27)  may  be  neg¬ 
lected  compared  to  those  varying  as  l/fi  ^  in  equation  (9.28). 

With  the  above  approximations  add  equation  (9.27)  and  (9.28)  to 
obtain  the  resultant  air  gap  vector  potential  due  to  the  rotor  current 


such  that 


V 


< 


/; 


■  j- 


;! 


^  4 

:  v-  v.-:.  ,  . :  ■  .  ■; .  <;/  ;  .  *  =.  . 

■ 


r 


, 

> 


113 


A^(T  9>t) 


M 


\1\  d 


0  1  C°S  t  +  %) 


-^o  (1  +  £.)  cose®  +%)  t 


-°tfr 


00 

'W  f^/f  c^ec°s[|^-f(i+2v;]e  sit/3„t 

U  ...9.29 


V—  0 


where 


4  — 


3  = 


Tb-4 


>J  _  6j5  L  Ij_ 


^l=/[(Rrf+Ct03Lr^] 


4  = 


 £v 


LtC 


The  resultant  vector  potential  due  to  the  stator  currents  only  is 
obtained  by  adding  equations  (9.3)  and  (9.4) 
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9.4  Transient  Field  Intensities 


Since  \o(v j  «/Sv  »  "the  electric  field  intensities  are 
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9. 5  Transient  Torque  on  the  Rotor 

From  section  (6.4)  the  torque  on  the  rotor  is 
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This  is  the  rate  at  which  electromagnetic  angular  momentum  is  flowing 
inward  through  the  rotor  surface 
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9 . 6  Storage  of  Electromgnetic  Angular  Momentum  in  the  Air  Gap 

From  equation  (6.12)  the  electromagnetic  angular  momentum  stored 
per  unit  axial  length  of  the  air  gap  is 
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Evaluating  the  integral  with  respect  to  £> 
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In  terms  of  the  above  constants, 
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This  is  the  electromagnetic  angular  momentum  stored  in  the  air  gap 
as  a  function  of  time. 

9 . 7  Stator- Rotor  Torque  Unbalance -and  Torque  on  the  Stator. 

From  equation  (6 . 16J,  Hie  Stator- Rotor  Torque  Unbalance  is  given  by 
Rate  of  change  of  stored  angular  momentum  in  the  air  gap. 
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The  above  equation  is  true -because-^  and  l°fyl  are  much  smaller 
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This  is  the  rate  at  which  electromagnetic  angular  momentum  is  trans¬ 
mitted  across  the -air  gap  surface  at  r  =  b • (that  is  stator  surface). 
9 . 8  Torque  Expressions  for  small  values  of  the  Time. 

(A)  Consider  the  rotor  torque  given  by  equation  (9.39).  Define 
the  Time  Constants  J  and  j  as 
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two  series  cancel  each  other  termwise.  Therefore, 
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This  follows  from  the  fact  -that  Tv  is  independent  of  v  . 
The  above  series  -may  -be  written  as 
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This  is  a  well  known  Fourier  series  whose  sum  is 


£?n=  I,  CK«f<7T 
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It  follows  that 
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( B )  Cons  ider  now  the  stator  torque  given  by  equation  (9.46).  For 
small  values  of t due  to  the- same -reasoning  as  for  the  rotor  torque, 
the  first  two  terms  and -last -two  series -may  be  omitted.  Furthermore, 
it  may  be  shown  that  the -magnitudes  of  the  third  and  the  fourth 
terms  are  negligible  in  comparison  -with  the  magnitude  of  any  other 
term.  Hence, 
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Now  the  series 
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This  is  again  a  Fourier  series  and  may  be  easily  summed  by  using  the 
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'Geometric  series  Method'  .  The  sum  is 
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9.9  Torque  Expressions  for  Large  Values  of  the  Time 

Since  <7  >>  7  ,  for  1  <t<J\  the  contributions  of  the  series  terms  in 
v5  v  5 

equations  (9.39)  and  (9.46)  are  negligible.  Also  for  reasons  stated  in  the 
last  section,  the  third  and  the  fourth  terms  in  equation  (9.46)  may  be 
omitted .  Hence , 
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Also  it  follows  that  for  t>>  7 
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=  steady  state  torque 
9.10  Numerical  Results 

(A)  The  Time  Constants:  From  equation  (9.23) 
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Since  Rx«fJt>ci  one  may  write 


Hence , 


.9.51 


Furthermore ,  using  equation  (9.22)  one  obtains 
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In  general,  therefore, 

T»TV 


The  time  Constants  are  calculated  as  functions  of  the  Air  gap 
length  for 

a  =  1  m  =  Rotor  Radius 

Rp  =  1.72  XI cA=  Rotor  surface  resistivity 
The  results  are  shown  as  graphs  in  Fig.  (92). 

(B)  Magnitude  of  the  Torque  components:  Let  equations  (9.47), 
(9.48),  and  (9.50)  be  written  as 
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...9.55 


Where 


Tn  =  2Fu)  La2b2/cd 
1  s 


...9.56 


T2  =  7Ta)sLb2  2cd 


9.57 


...9.58 


For  the  values  of  a  and  given  in  part  (A)  of  this  section, 

the  magnitudes  of  T-^,T2  and  (these  are  the  torque  magnitudes 


s  * 

normalized  with  respect  to  1^)  are  calculated  as  functions  of  the  air 


gap  length.  The  results  are  plotted  in  Fig.  (9.3).  Note  that  T-^  gives 
the  normalized  value  of  the  rotor  torque  at  t  =  0+. 

(C)  Normalized  Torque  Versus  Time  (t  small):  For  the  above  values 


the  periodic  variation  Tr  (t)  and  Ts  (t)  given  by 
^  em  em 


of  a  and  R 


equations  (9.53)  and  (9.54)  respectively  are  shown  in  normalized  form 
in  Fig.  (9.4).  These  curves  may  be  used  to  determine  the  differential 
torque  impulse  on  the  stator  and  rotor  (i.e.  the  difference  in  torque 
integrated  over  the  time  of  the  transients). 

9 .11  Discussion 

(A)  Discussion  on  the  Solution:  To  advance  the  concept  of  electro¬ 
magnetic  angular  momentum  flow  in  an  induction  motor  a  transient 
problem  has  been  investigated.  The  solution  of  the  field  equations  is 
obtained  in  the  exact  sense.  In  the  equations  giving  the  electromag- 
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g  -  Air  Gap  Length  in  m. 

Fig.  9.2  -  Time  Constants  Versus  Air  Gap  Length 
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netic  quantities  only  those  terms  are  omitted  whose  magnitudes  are  small 
compared  to  the  magnitudes  of  the  terms  retained.  It  may  be  shown  that 
under  the  quasi-static  approximation  the  series  terms  will  be  missing 
from  all  these  equations  and  such  a  solution  will  not  satisfy  the  initial 
conditions  of  rotor  vector  potential  waves.  As  a  result  the  initial 
values  of  the  torque  will  not  be  properly  treated. 

(B)  Discussion  on  Amentum  Effects:  Consider  the  expressions 
for  the  rotor  torque  given  by  equation  (9.39),  (9.47)  and  (9.49).  These 
equations  may  be  interpreted  as  follows: 

From  t  =  0  the  rotor  starts  absorbing  electromagnetic  angular 
momentum  stored  in  the  air  gap  by  the  stator  field.  The  combined  effects 
of  all  the  terms  (including  those  neglected)  in  equation  (9.39)  should 
give  a  continuous  rate  of  change  of  absorption  (this  removes  the  dis¬ 
continuity  of  T^m(t)  in  equation  (9.47)  at  t  =  0)  such  that  at  t  =  0+ 
the  rate  at  which  the  momentum  is  crossing  the  rotor  surface  inward  is 
approximately  given  by  equation  (9.47).  As  this  function  decreases 
the  effect  of  other  terms  increases.  At  t  =  2g/c  +  0++  (Fig.  (9.4)) 
the  rotor  starts  giving  out  electromagnetic  momentum.  The  process  of 
absorption  and  rejection  continues  in  cyclic  order  and  for  t  >7vthe  rate 
is  approximately  given  by  equation  (9.39),  which  would  be  obtained  under 
static  approximation. 

Consider  now  equations  (9.46)  and  (9.48),  which  are  the  expres¬ 
sions  for  the  stator  torque.  A  close  examination  of  these  equations 
using  typical  numerical  values  will  show  that  for  very  small  values 
of  t  (of  the  order  of  fraction  of  a  y^s)  the  stator  torque  is  approx¬ 
imately  zero.  This  is  an  indication  that  the  torque  on  the  stator  comes 
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from  a  delayed  action,  which  is  in  support  of  the  momentum  concept.  At 
t  >  the  stator  and  rotor  torques  are  equal  in  magnitude  and  are  given 
by  equation  (9.49). 

From  the  above  discussions  it  follows  that  the  quasi-static 
approach  can  not  give  proper  account  of  the  torque  at  the  first 
instants  (t  <^)  during  the  course  of  the  transient.  The  time  functions 
associated  with  the  large  roots  of  the  wave  system  are  responsible  for  the 
initial  flow  of  the  momentum. 

(c)  Discussion  on  the  Magnitude  of  the  Torque  Components:  The 
minimum  air  gap  length ^  of  an  induction  machine  with  a  =  1  m  is  about 
.0017m.  From  Fig.  (9.3)  it  follows  that  T-^/Tg  giving  the  transient 
torque  at  t  =  0+  to  the  steady  state  torque  for  Rp  =  1.72  ( 10“ 14 )  is  about 
1.3.  This  ratio  increases  if  ^decreases  (Tq  constant,  T3  decreases  with 
Rr  almost  in  direct  proportion) .  In  a  real  machine  normally  Rr  would  be 
larger  and  Tq/Tg  would  be  smaller  than  the  above  value. 
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CHAPTER  10 


CONCLUSIONS 

The  law  of  conservation  of  momentum  has  been  interpreted  for  an 
electromagnetic  system.  This  has  led  to  a  concept  of  a  conversion  of 
electromagnetic  angular  momentum  to  mechanical  angular  mometum  to 
produce  torque.  With  a  view  to  apply  this  in  the  theory  of  electrical 
machines,  an  expression  has  been  obtained  in  terms  of  Maxwell’s  stress 
tensor  and  electromagnetic  angular  momentum  density  to  calculate  the 
electromagnetic  torque  on  a  body  in  free  space. 

The  electromagnetic  field  problem  of  double  cylindrical  machines 
has  been  reduced  to  finding  the  transient  and  the  steady  state  solutions 
of  a  vector  wave  equation.  A  general  solution  of  this  equation  has 
been  developed  by  using  the  method  of  the  Laplace  transform.  This 
exact  solution  may  be  utilized  for  the  harmonic  (space  and  time) 
analysis  of  certain  types  of  double  cylindrical  machines. 

The  machine  volt-ampere  equations  have  been  obtained  by  inte¬ 
grating  the  electric  field  intensities  and  the  torque  equation  has  been 
deduced  from  the  momentum  concept.  The  quasi-static  field  quantities 
and  the  resulting  machine  equations  of  motion  have  also  been  derived 
from  the  low  argument  approximations  of  the  Bessel  functions  in¬ 
volved  in  the  exact  equations.  These  results  are  in  agreement  with 
those  found  by  other  methods. 

A  detailed  wave  theory  analysis  of  an  induction  machine  has  been 
developed  in  the  exact  sense  to  study  the  production  of  torque  in 
terms  of  a  transfer  of  electromagnetic  angular  momentum  across  the 
air  gap.  It  has  been  shown  that  the  quasi-static  approximation  gives 
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completely  satisfactory  quantitative  results  (accuracy  -  one  part  in 
about  10l0)  in  the  case  of  the  machine  under  steady  state  conditions. 
However,  in  the  transient  case  such  an  approximation  cannot  provide 
information  about  the  torque  for  the  first  instants  (of  the  order 
of  yU  seconds)  after  the  transient  has  started.  An  exact  solution  is 
necessary,  and  in  such  a  solution  the  time  functions  associated  with 
the  large  roots  of  the  electromagnetic  wave  system,  obtained  in  the 
form  of  infinite  series,  are  responsible  for  the  initial  flow  of 
electromagnetic  angular  momentum. 

No  doubt,  one  is  fully  entitled  to  explain  the  production  of 
torque  in  the  induction  machine  from  the  concept  of  a  flow  of 
electromagnetic  angular  momentum  and  its  conversion  into  mechanical 
angular  momentum. 
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APPENDIX  1 


In  this  appendix  it  will  be  shown  that 

£0|!(TT£J  -E x(vxEj] 


+^oi$  iy-  b )  -Hx^xtyj 


=  M,fP  =  AH5l4  V"*8,  <  =  1,2,3  ...Al.l 


up 


where 
Proof : 


..A1.2 


It  is  known  that 


ana 


=  £, 


E0  [£(  V-E)  -  E  xOxej] 

[  £  ( 7-E; -f  (  E  v)-t]  -  i  [E.(Ef)  +  £.(?vJ i] 


. . .Al. 3 
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. . .A1.4 


Also  it  can  be  shown  that 
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Furthermore 


CCu-i;) U]-  V  =  V  (u.w)  =  (u^-)V 


»  . .  A1  o  7 


Using  the  facts  given  by  the  last  three  equations,  if  equa¬ 
tions  (A1.3)  and  (A1.4)  are  added,  the  desired  result  will  follow. 
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APPENDIX  2 


In  this  appendix  Maxwell's  equations  will  be  solved  in  terms 
of  suitably  defined  potential  functions. 

Since  V . §  =  0 ,  one  may  write 

5  =  V  x A  ...  A2.1 

v 

From  VxE  =-  —  and  equation  (A2.1)  it  follows  that 


8t 

f  =_  _ 


v$ 


at'1  ...  A2. 2 

where  A  and  0  and  electromagnetic  vector  and  scalar  potentials  respect¬ 
ively  . 

Substituting  equation  (A2.1)  and  (A2.2)  in 

v-a-Vtr 

V.S  =  f> 

and  using  the  fact  that  for  homogeneous  and  isotropic  medium  sta¬ 
tionary  with  respect  to  the  reference  system  56  =  €-56  and  5  =  yf5, 
one  obtains 


VxVx£  +  uev|l  +]szd—^  =  yjp 

at2  f 


...  A2.3 


V2$  +  V.  =-  P/e  •••  A2*4 

For  identical  boundary  conditions  all  particular  solutions  of 
(A2.3)  and  (A2.4)  would  lead  to  the  same  electromagnetic  fields.  Since 
A  and  $  are  not  uniquely  determined  by  (A2.1)  and  (A2.2.),  a  further 
restriction  may  be  applied  in  order  to  simplify  the  equations.  One  such 
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restriction  is 


V 


1 


+ 


V  E 


9_$ 
3  t 


=  0 


. . .  A2 . 5 


which  is  known  as  the  Lorentz  condition. 

Substitution  of  equation  (A2.5)  in  equations  (A2.3)  and  (A2.4) 
yields 


y  e 


32A 

3t2 


V2$ 


-  ye 


9  2$ 

9t2 


=  -  P/e 


...  A2 .6 


. ..  A2.7 


For  the  class  of  problems  dealt  with  in  this  thesis,  p=  0 
in  all  free  space  regions  and,  then  it  is  possible  to  choose  $  =  0. 


This  is  equavelent  to  the  condition 

V.A  =  0 
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APPENDIX  3 

The  evaluation  of  the  residues  in  chapter  9  involves  the  use 

of  the  roots  of  the  function 

[K1  (kb)  I  1(ka)  -  I,1(kb)K  (ka)] 

Rv>  +  Q  [K1  (ka)  I'  (kb)  -  I'  (ka)K'  (kb)]  =  0 


where 


...  A3 .1 


and, 


=  Rotor  surface  resistivity,  k  =  — 


Q  =  y  c  =  Constant  such  that  Q  >>  R 

v  o  ■£ 


The  required  roots  will  be  determined  as  follows : 

50 

(a)  Large  Roots:  Normally  the  large  roots  of  the  Bessel 
functions  and  certain  other  related  functions  are  obtained  by  using 
the  asymptotic  expansions  of  these  functions  valid  for  large  argu¬ 
ments  in  the  required  sector  of  the  complex  plane.  This  principle 

51 

and  the  properties  of  asymptotic  series  may  be  applied  to  find 

the  large  roots  of  equation  (A3.1).  The  modified  Bessel  functions 

50 

of  order  one  have  the  asymptotic  expansion  , 
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where , 


A1  *2 

U(Z)  =  1  -  —  +  -4 

z  l 
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A3  A4 

T  "4 

z  z 


A5  + 

T  “• 

z 

...  A3. 4 


where, 


A.  A0  A0  A.  A,- 

V(Z)  =  i  +  ^+  4  +  4  +  4  +  4  + 

z  z  z  z 


...  A3 . 5 


Let 
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A  14175  a  _  14175.(77) 
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...  A3 . 6 


Differentiating  equation  (A3. 6),  one  obtains 
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with 


B  =  A  + 
n  n 


(2n-l) 

2 


Similarly 


and 


K* 1(z) 


-z+0 ' 
e 


Let 


n  >  2 


. . .  A3 . 10 


•  •  • 


A3. 11 


|  =  Pi  ka=z 
/.  kb  =  Pz 

Using  these  relations  with  equations  (A3. 6),  (A3. 7),  (A3. 10)  and 
(A3. 11),  it  may  be  shown  that 


K,1(kb)I1(ka)  -  If1(kb)K1(ka) 


-  — —  cosh  \(P- l)z  +  9l-e  1  (p) 1 
2  /pz 


...  A3. 12 


K,1(ka)I'1(kb)  -  I ,1(ka)K’1(kb) 


1 


sinh  r(P-l)2+  e*  -  0*  (p)l  ...  A3. 13 

2/7%  1  1  1  J 

~0  *  (P)  +0  ’  (P) 

where,  e  1  rv  and  e  1  r  are  given  by  equation  (A3. 8)  and 


(A3. 9)  respectively  with  z  replaced  by  pz. 
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Substitution  of  (A3. 12)  and  (A3. 13)  into  (A3.1)  reduces  it  to 


cosh  [CP  — l)z'  +  0  -  0'  (P)] 

sinh  [(p-l)z  +  8^  -  6 1  x  (pi) 


where 

R  =  R^/Q 

...  A3 . 15 

Let 

(p-l)z  =  mz 

and 

*  =  0-l  -  01!  <P) 

»  =  S'  -  0'  (P) 

Hence , 

cosh(mz+t)  _ 

"  —  —  I\ 

sinh(mz+$) 

from  which 

...  A3 . 16 
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B  B  ,  B  0 

c  =  -a  -  R  +  J2i|  B,  +  . 
n  />n  «n-l  1  aFl-2  2 


pn  pl 


P1 


.  +  (-l)nB 


n 


. . .  A3 . 18 


B1 
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...  A3. 19 
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.  A3. 20 


Substitution  of  these  relations  in  equation  (A3. 16)  leads  to 
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Equation  (A3. 21)  may  be  expressed  as 
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Such  that 


a 

F  =  t-2. 
o  b 
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Solution  of  Equation  (A3 .22) 

Approximate  solution  of  equation  (A3. 22)  will  now  be  obtained  by 

numerical  method. 

First  order  approximation: 

Let  e  =  F 
o 

R  —  1 

The  quantity  Fq  =  ^  +  ^  is  a  negative  number  with  magnitude  being 
slightly  less  than  1.  Hence 

1  -  R 


2mz  „ .  /  „ 
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1  +  R 


Solving  for  z,  one  obtains 
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Second  order  approximation 
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Re(2)  =  i-  m  (W  + 


lv  2' 


...  A3. 26 


Im(z)  =  6  + 
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Proceeding  in  the  same  way  it  may  be  obtained  that 


Re(z)  =  —  In 
4m 


M  M 

(M  +  -i  +  it  +  )2 
°  62  6 4 


M  M  M 
1  3  6  7 

+  (—  +  —  +  —  + 

6  8 3  85  ' 


...  A3. 30 


Im(z)  =  3  + 


h  N3  N5 
8  +  B3  8 5 


...  A  3.31 


.> 


as.iA  . 


i 


h- 


■ 


J.- 


■ 


- 1 


fc  -  t 7 


k, 


' 


- 

<* 


X 

•  -V  s' 


--  • 


. 


-  ^  . 

•  ,  i 

, 


OC.£A  ... 


/ 


pH  .  ^  ■ 

' 


( 


■ 


/" 


.  -V 


. 


148 


where 


\ 


Mo  =  Fo>  M2  ='F2’  M4  =  (2F2A  +  V 


Mi  =  Fr  M3  =  -  (FlA  +  F3) 


M,-  =  (F,A  +  F,B.  +  3F0A  +  Fc) 

O  1  11  0  o 


Nl  =  5--  N,  =  -  ^  +  F1A  +  F3n. 
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Equations  (A3. 30)  and  (A3. 31)  representing  the  real  and  imag¬ 
inary  parts  of  z  will  be  sufficiently  convergent  to  give  the  large 
roots  quite  accurately.  The  earlier  roots  will  be  approximate. 

Since  ka  =  =  z,  the  large  roots  of  equation  (A3.1)  are. 


therefore,  given  by 
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in  conjugate  pairs. 

(B)  Samll  Roots:  The  small  roots  of  equation  (A3.1)  may  be 

obtained  by  using  the  ascending  series  for  the  Bessel  functions 

48 

involved.  These  series  are 


Il(z)  2  +  16  +  **• 


r=o 


7^Z)1+2r 

r!  r(r+2) 

I 


...  A3 . 35 


K,  (z)  =  —  +  i-z'lnz  +  . . 
1  z  2 


T  •  I  (2)  -  I  (Z) 
Lim  -n  n 


n  ->  1 


sin  nil 


...  A3 . 36 


Using  only  the  1st  term  of  each  of  the  above  series  as  in 
section  (7.2),  it  may  be  shown  that 


Lim  [Kf1(kb)I1(ka)  -  I'  (fcb^Cka)] 

[K!  ,(ka)r:(kh)  -  I’1(ka)K'1(kb)J 


s  o 


sa  [b/a  +  a/bl 
c  [b/a  -  a/b] 


. . .  A3 . 37 


Let 


T  =  „  =  Lb/a  +  a/b] 
r  y  o  '  j^b/a  -  a/bj 


. . .  A3. 38 


Therefore,  the  smallest  zero  of  the  left  hand  side  of  (A3.1)  is 

s  =  -  R  /L  ...  A3 . 39 

r  r 


Using  more  terms  of  the  series  in  equations  (A3. 35)  and  (A3. 36) 
estimation  shows  that  for  —  =  P,  o{  P  _<  1,  the  absolute  magnitude  of 

cl 

the  next  higher  zero  of  equation  (A3.1)  is  greater  than  c  (magni¬ 
tude  of  the  velocity  of  light).  If  S£  is  this  zero, 
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k 


>  1 


Since  for  this  condition  the  asymptotic  series  shown  in  equa¬ 
tions  (A3. 2)  through  (A3. 5)  are  valid,  the  numerical  value  of  may 
be  calculated  as  in  section  (A)  of  this  appendix,  and  is  approximately 
equal  to 

[l  -  R  hi  c| 

c  ,  L  r'  o  J  .  lie 

S2  "  2(b-a)  [1  +  Rh c]  3  2(b-a)  ...  A3. 40 

with  another  zero  occurring  in  conjugate  pair. 

This  completes  the  determination  of  all  the  roots  of  equation 
(A3.1),  with  their  first  order  approximate  values  given  by  equations 


(A3. 34)  and  (A3. 39). 
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APPENDIX  4 

In  this  appendix  it  will  be  proved  that  when  a<  r  ^  b • and 
t  >  0 ,  the  integrals 


h  = 


tt  rr,  x  st  , 
co  X,:  (r,s)  e  ds 
si. 


s(.sZ  +  «_2)  [Rr  +  uoc  X^s)]  . . .  A4.1 


o  1 


h  = 


tt  rr,  v  st  , 
r  X^  (r,s)e  ds 


(s^  +  w  VR  +  y  cX^  (s)l  ...  A4.2 

s  i  L  r  o  1  J 


vanish  in  the  limit  over  the  arc  FGE  of  Fig.  (9.1). 

Proof :  For  s  tending  to  infinity  and  co  small 

s 


,  2  ,  2  ■  3 

s(s  +  CO  )  =  s 
s 


and 


2  2*2 
s  +  co  =  s 
s 


Furthermore,  since  R  =  0  and  yu  c  >>  o  for  a  path  where 


X-^Cs) 


>  o 


the  function 


+  y  qc  X^^Cs)  =  /^C  X^r(s) 


The  problem  is  thus  reduced  to  show  that  for  a  ^  r  b  and  t  >  o , 
the  integral 


I  = 


r  X^Crs) 

r  1  1.  5  st  J 

- e  ds 

rr,  , 

y  cs:  X^  ( s ) 


...  A4. 3 


m  =  2  or  3 
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vanishes  in  the  limit  over  the  required  arc  such  that  on  this  path  , 

Xn  SS,  v  >  0. 

1  (s) 

Using  the  asymptotic  expansions  of  the  Bessel  functions  and 
following  the  procedure  as  in  appendix  3 ,  it  may  be  shown  that  for  s 
tending  to  infinity. 


X^^OjS)  = 


[K'-^Csb/c)  1^  (sr/c)  -  If^  (sb/c)  K^(sr/c)^) 
[KUCsa/c)  IU  (sb/c)  -  I ' (sa’  /c)K' ,  (sa/c)] 


=  /—  cosh  s(b-r)/c 
siTih  s(b-a)/c 


...  A4.4 


and 


_  [K'^Csb/c)  1^  (sa/c)  -  I’-^Csb/c)  K^(sa/c)] 

^1  ^ s ^  ^KT^(sa/c)  IT-^(sb/c)  -  I'-^Csa/c)  K'^(sb/c)J 


[cosh  s(b-a)/c] 
[sinh  s(b-a)/cj 


...  A4.5 


From  which 


X^^CrjS) 

X^Cs) 


/a-  cosh  [s(b-r)/c] 
/  r  cosh  [s(b-a)/c] 


Choose  Hie  radius  of  the  circular  arc  R  = 


nvc 

b-a 


...  A4 . 6 


with  points  on 


the  circumference  given  by  s  =  ReJ  ,  -  II  _<  <p  _<  n.  From  equations  (A4.5) 


and  (A4.6)  only  a  few  steps  will  be  needed  to  prove  that  over  the 

TV, 


contour  defined  by  R,  the  function 


X^  ( s ) 


>  o  and  the  factor 


X1rr(r,s) 


X^s) 


cosh  [s(b-r)/c] 

cosh  [s(b-a)/cj 


<  c 
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showing  that  it  is  bounded.  Furthermore,  the  location  of  the  vth  zero 
of  the  function  R  i  y  qcx^  (s)  =  0,  defined  by  equation  (A4.34)  is  at 
a  distance  of 


R 

r _ 

y o(b-a) 


+ 


1  2 

TpC  (  2V-1)  TtVC 

2(b-a)  <  b-a 


from  the  origin.  Hence,  the  arc  does  pass  through  this  zero  and  all 
other  zeroes  having  smaller  magnitudes  ( including  that  at  s  =  -  r  )  . 


L 


Therefore,  the  absolute  value  of  I 


I 


r  v  rr,  s  st 

rxx  (r,s)e 


^cs^Cs) 


ds 


<  R  R 
m 


-m+1 


r 


etRcosH+ 


where,  the  relations 

1 

=  R  m  1  ds  1  =  Rd  t  and 

st 

e 

m 

s 

;  1  1 

^tRcoscf) 


have  been  used  and  R  is  a  constant  which  is  bounded. 

m 


Now  for  the  path  FF1  and  EF’ 

tRcos^qt;  ,y*>o 


so  that 


IFF’ 


IEE ' 


V2 

<  R  R“m+1ert  f  d4> 
m 


<}> 


1 


■a  ^-m+i  n  .  -l  T 
=  R  R  e  sin  jr 
m  R 


rb-ai^-l  Tt  .  -1  Y'(b-a) 
- — J  e  sin  - t— — 

m  L  it  vcJ  tt  Vc 


...  A4.7 


V 


1  ... 


■  ... 


V--'  if. 


.... 
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-v^  ? 
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V-  -V 
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Hence 


Lim 


IFF' 


=  Lim 


^EE’ 


=  o 


v  ->■  00 


V  oo 


...  A4 . 8 


For  the  path  F'G  and  EfG 


IF’G 


ZE’G 


R  R 
m 


-m+1 


R  R 
m 


e  <34, 


77/2 

ttV2 


-m+1 


e  -tRsin*  ^ 


R  R 
m 


o 

TT  /  2 


-m+1 


r 


e  -2tR*/1'  <34 


U 


-r,  -n-Hl  TT-  . 

<  RmR  2t  >  1  "  ° 


=  R 


r  b-a 


m  Lfj-vc  -*  2t 


m 

]  W’  *  >  ° 


.  A4.9 


This  shows  that 


Lim 
v  -* 


Lim 
v  -*  °° 


. ..  A4.10 


Therefore,  as  the  integrals  of  equations  (A4.1)  and 


(A4.2)  vanish  over  the  arc  FGE  of  Fig.  (9.1). 
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